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NUMERICAL  AND  ALGEBRAIC 


Gain  in  decibels  of  P2  relative  to  Pi 

G  =  101og10(P2/Pi). 


To  within  two  percent 

(2tt)1/2  «  2.5;  *2  «  10;  e3  «  20;  210  a  103. 


Euler-Mascheroni  constant1  7  =  0.57722 

Gamma  Function  T(s  +  1)  —  xT(x): 

T(l/6)  =  5.5663  r(3/5)  =  1.4892 

T(l/5)  =  4.5908  T(2/3)  =  1.3541 

T(l/4)  =  3.6256  T(3/4)  =  1.2254 

T(l/3)  =  2.6789  T(4/5)  =  1.1642 

T(2/5)  =  2.2182  T(5/6)  =  1.1288 

T(l/2)  =  1.7725  =  v^F  ^(1)  =  1.0 


Binomial  Theorem  (good  for  |x|  <  1  or  a  —  positive  integer): 

OO 

„  v — /a\  jl  a(a  —  1)  2  <*(“  ~  1)(Q:  —  2)  3 

(i +x)  (*) 1  = 1 + as + “4H1  +  ~ — 1  +*-- 


Jc=0 


Rothe-Hageu  identity2  (good  for  all  complex  z,  y,  z  except  when  singular): 


Ex  tx  +  kz\  y  /y  i-  it 1  -  k;z\ 

x  +  kz'  k  '  y  +  (n  —  k)z  '  n  —  k  ' 


y  +  (n  -  fc)z> 


fc=o 


x  +  y  fx  +  y  +  nz\ 
x  +  y  +  nz  '  n  ' 


Newberger’s  summation  formula3  [good  for  n  nonintegral,  Re(a  +  /3)  >  -1]: 

V  Lt  _  _JL -Ja+^(Z)J0_^(Z). 

/  j  n  -f  y.  sin  n* 
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VECTOR  IDENTITIES4 


Notation:  /,  g ,  are  scalars:  A.  B,  etc.,  are  vectors:  T  is  a  tensor:  /  is  the 
unit  dyad. 

(1)  A'BxC  =  AxB'C  =  B'CxA  =  BxCA  =  CAxB  =  CxAB 

(2)  A  x  (B  x  C)  =  (C  x  B)  x  A  =  (A  •  C)B  -  (A  •  B)C 

(3)  A  x  (B  x  C)  +  B  x  (C  x  A)  +  C  x  (A  x  B)  =  0 

(4)  (A  x  B)  •  (C  x  D)  =  (A  •  C)(B  •  D)  -  (A  •  D)(B  •  C) 

(5)  (A  x  B)  x  (C  x  D)  =  (A  x  B  •  D)C  -  (A  x  B  •  C)D 

(6)  V(/5)  =  V(gf)  =  fVg  +  gVf 

(7)  V  •  (/A)  =  fV  ■  A  +  A  •  V/ 

(8)  V  x  (/A)  =  /V  x  A  +  V/  x  A 

(9) V.(AxB)  =  B-VxA  — A-VxB 

(10)  V  x  (A  x  B)  =  A( V  ■  B)  -  B(V  •  A)  +  (B  •  V)A  -  (A  •  V)B 

(11)  A  x  (V  x  B)  =  (VB)  •  A  -  (A  •  V)B 

(12)  V(A  •  B)  =  A  x  (V  x  B)  +  B  x  (V  x  A)  +  (A  •  V)B  +  (B  •  V)A 

(13)  V2/  =  V  ■  V/ 

(14)  V2A  =  V(V  •  A)  -  V  x  V  x  A 

(15)  V  x  V/  =  0 

(16)  V-Vx  A  =  0 

If  ei ,  e2,  es  are  orthonormal  unit  vectors,  a  second-order  tensor  T  can 
be  written  in  the  dyadic  form 

(17)  Ta'£(,iTv“i 

In  cartesian  coordinates  the  divergence  of  a  tensor  is  a  vector  with  com¬ 
ponents 

(18)  (V-T)i  =  Y^^dTji/dzj) 

[This  definition  is  required  for  consistency  with  Eq.  (29)].  In  general 

(19)  V  •  (AB)  =  (V  •  A)B  +  (A  •  V)B 

(20)  V  •  (/T)  =  V/-T+/V-T 
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Let  r  =  ix  +  jy  +  kz  be  the  radius  vector  of  magnitude  r,  from  the  origin 
to  the  point  x,y,  z.  Then 

(21)  V  •  r  =  3 

(22)  V  x  r  =  0 

(23)  Vr  =  r/r 

(24)  V(l/r)  =  -r/r3 

(25)  V  *  (r/r3)  =  4x*(r) 

(26)  Vr  =  / 


If  V  is  a  volume  enclosed  by  a  surface  S  and  dS 
unit  normal  outward  from  V, 


txdS  jrhere  n  Is  !&i 


(27)  L  dWf  =  J  dSf 

J 

-I, 


(28)  j  dVV  •  A  =  /  dS-A 
v  ./  s 


(29)  j  dW-T  =  /  dS-T 
v  J  s 


(30)  j  dW 
v 


x  A  =  J  dS  x 


(31)  /  dV(/V2g  -  yV2/)  =  /  dS  •  (/Vy  -  y V/) 


(32)  /  dV(A  •  V  x  V  x  B  -  B  -  V  x  V  x  A) 

=  /  dS  •  (B  x  V  xA  —  AxV  xB' 


If  5  is  an  open  surface  bounded  by  the  contour  C.  of  which  the  line 
element  is  dl. 


(33) 


dS  x  V/ 


dl/ 
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<W) 

(35) 

(36) 


L 

I, 


dS-V  X  A  = 

(dS  x  V)  x  A  =  < 

J 

dS-(Vf  x  Vg)  = 


d  1  •  A 


x  A 


fdg  = 


gdf 


DIFFERENTIAL  OPERATORS  IN 
CURVILINEAR  COORDINATES5 


Cylindrical  Coordinates 
Divergence 


V  •  A 


1 

r 


+ 


5a 


Gradient 


(V/). 


<9/. 

dr' 


(V/)*  = 


15/. 

r  d<i> 


(V/), 


Curl 


(V  x  A)r 


i  axr 

r  c?</> 


d  A* 
dz 


(V  x  A)d  = 


dAr 

dz 


dA; 

dr 


d± 

dz 


(V  x  A) 


1  d  1  dAr 

*  -  ~  —(rA<t>) - 


dr 


r  d(p 


Laplacian 


V2/ 


r  dr  \  dr  J  r2  d<P 2 


<9z2 


6 


Laplacian  of  a  vector 
(V2A),.  = 

r2  dp  r2 


(V2A)0  =  v2aC)  +  — 


2  dA, 


t2  dp 


( V2  A)r  =  V2Ax 
Components  of  (A  •  V)B 


dBr 

A*  dBr 

dr 

f  r 

dB* 

dr 

r  <9# 

dBz 

a*  asx 

dBr  A*B0 
dz  r 


dB *  AaBr 

+  At  — —  +  - 

dz  r 


dr  r  dp 


Divergence  of  a  tensor 

(V  •  T)r  =  ~~(rTrr)  +  ^  -  — 

ror  r  oz  r 


(V  •  T),  =  i  i-(rr,,)  +  i  ^ 

r  Or  r  dp  dz  r 


, _  i  a  .  i  arOJ  at., 

(V-  T),  =  —  ~(rTP, )  +  -r“ 

r  or  r  dp  dz 
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Spherical  Coordinates 

Divergence 


V  •  A  = 


Id,  Id 

~  —  (r'Ar)  +  — -  —  (sin#A*) 
r-  dr  r  sm  #  do 


1  0Ao 
r  sin  8  d( ? 


Gradient 

(V/)r 


£/. 

dr  ’ 


(V/),  = 


1  £/. 

7  d#  ’ 


(v/)0  = 


i  a/ 

r  siu  5  d<£ 


Curl 


(V  x  A)r 


d 


r  sin  8  88 


-—(sin#A*)  - 


1  8Ae 
r  sin  8  dp 


(V  x  A)e  = 


1  8Ar 
rsin#  dp 


id_ 
r  8r 


(rA*) 


.  1  d  ,  .  .  1  8 Ar 

(V  x  A)*  =  --fc,(TA.)  -~-gf 


Laplacian 

V2/ 


_l__d_  /  2  d/\  1  8  f  d/\  1  d2/ 

r2  dr  \  dr  J  r2  sin  8  89  \  88  )  +  r2  sin2  9  dtp2 


Laplacian  of  a  vector 


(  V2  A)r  =  V2/4r 


2  Ar 
r2 


2  dA$  2cot#A?  2  dAd 

r2  d#  r2  r2  sin  #  dp 


(V  A)#  =  VM,  +  — 


2  8 A, 


Aj  2  cos#  dAd 


r2  88  r2  sin2  9  r2  sin2  9  8<t> 


,  ,  A*  2  dAr  2  cos  8  8Aq 

(V2A)*  =  V2A*  -  - -Vr  +  -T— r—  + 


r*  sin 


in2#  r2  sin  #  dp  r2  sin2  #  d<£ 
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Components  of  (A  •  V)B 


(A  •  VB)r 


Ag  dBr 
r  d8 


A0  dBr  AgBg  -r  A0B0 

r  sin  8  d<p  r 


,  A  _  4  dBg  t  Ag  dBg 

(A.VB),=^— +  -  — 


A±  dBg  AgBr 

4 - — - -  4 - - — 

r  sin  8  dtp  r 


cot  8AoB0 


r 


(A  •  VB)*, 


^  dB &  Ag  dBe j _ A, 

dr  r  d8  r  sin  8  dtp 


dBj,  t  A^Br  ,  cot  8Ac,Bg 

“r  '  ~T 


Divergence  of  a  tensor 


(V  •  T)r 


i  a  .  .  .  ,  i  d 

-9  r  +  r  sin  9  d6 


r2  dr 


(sin  8Tgr) 


1  dTtr  Tgg  ■+■  T<j,& 

rsin0  dtp  r 


^■T)>  =  ^j-r(r^)  + 


d 


r  sin  9  d9 


—  {sin  9Tgg) 


1  dT^g  Tg  r 

r  sin  9  dtp  r 


cot  9Ta<i> 


r 


(V-  T)*  =  ^4-(r2T^)  + 
r2  dr 


r  sin  9  d9 


— —  (sin  9Tgt) 


1  dTot> 
r  sin  9  dtp 


Tor  ,  cot  9T<i>g 
4"  T 

r  r 
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DIMENSIONS  AND  UNITS 


To  get  the  value  of  a  quantity  in  Gaussian  units,  multiply  the  value  ex¬ 
pressed  in  SI  units  by  the  conversion  factor.  Multiples  of  3  in  the  conversion 
factors  result  from  approximating  the  speed  of  light  c  =  2.9979  x  1010  cm/sec 
«3  x  1010  cm/sec. 


Dimensions 

Physical 

Sym- 

SI 

Conversion 

Gaussian 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

Capacitance 

C 

9  9 

t2q- 
ml 2 

l 

farad 

9  X  1011 

cm 

ml'2l3'2 

Charge 

9 

9 

t 

coulomb 

3  x  109 

statcoulomb 

_ 1/2 

Charge 

9 

m  ' 

coulomb 

3  x  103 

statcoulomb 

P 

l3 

i3r-t 

density 

/m3 

/cm3 

Conductance 

tq 2 

i 

siemens 

9  X  1011 

cm/sec 

mi2 

t 

Conductivity 

a 

tq2 

ml3 

l 

t 

siemens 

/m 

9  x  109 

sec-1 

Current 

q 

ml'2l3'2 

ampere 

3  x  109 

statampere 

*  t  * 

t 

t2 

1/2 

Current 

density 

9 

m  ' 

ampere 

/m2 

3  x  10s 

statampere 

/cm' 

l2t 

l'/2t2 

Density 

P 

m 

l3 

m 

l3 

kg/m3 

10-3 

g/cm3 

Displacement 

n 

9 

m1'2 

coulomb 

12‘x  x  105 

statcoulomb 

l2 

l'/2t 

/m2 

/cm2 

Electric  field 

F 

ml 

rn1'2 

volt/ m 

-  x  10"4 

3 

statvolt/ cm 

t2q 

l'/2t 

Electro- 

£, 

ml2 

ml/2  [1/2 

volt 

i  x  10-2 

statvolt 

t2  q 

motance 

Emf 

t 

3 

ml2 

ml2 

Energy 

U,W 

t 2 

t2 

joule 

10  ‘ 

erg 

Energy 

density 

w,  e 

m 

It2 

m 

It2 

joule/m3 

10 

erg/cm3 

10 


Dimensions 

Physical 

Svm- 

SI 

Conversion 

Gaussian 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

ml 

ml 

Force 

F 

— 

— 

newton 

10° 

dyne 

t- 

t 2 

Frequency 

/>  v 

1 

t 

1 

t 

hertz 

I 

hertz 

ml 2 

t 

1  -li 

-  x  10 

9 

Impedance 

Z 

tq2 

7 

ohm 

sec/cm 

ml2 

t 2 

Inductance 

L 

henry 

-  x  10-'1 

se  c'/cm 

q2 

i 

9 

Length 

l 

i 

i 

meter  (m) 

102 

centimeter 

(cm) 

Magnetic 

H 

q 

mx'2 

ampere- 

4t  x  10“3 

oersted 

it 

intensity 

turn/m 

Magnetic  flux 

ml2 

ml'2l3'2 

<J> 

B 

tq 

m 

t 

_ ,1/2 
m  ' 

weber 

tesla 

10s 

104 

maxwell 

gauss 

Magnetic 

tq 

2^/2  £ 

induction 

Magnetic 

l2q 

m1/2 15/2 

ampere-m2 

103 

oersted- 

m,  fi 

t 

t 

moment 

cm3 

Magnetization 

M 

q 

„ ,1/2 
m 

ampere- 

turn/m 

10~3 

oersted 

it 

Magneto- 

M , 

q 

ml/2  jl/2 

4t 

gilbert 

ampere- 

t 

10 

motance 

Mmf 

t2 

turn 

Mass 

m,  M 

m 

m 

kilogram 

103 

gram  (g) 

(kg) 

Momentum 

ml 

ml 

103 

g-cm/sec 

p.P 

t 

t 

kg-m/s 

Momentum 

density 

m 

Wt 

m 

Wt 

kg/ m2-s 

10_1 

g/cm2-sec 

ml 

1 

Permeability 

A* 

1 

heilry/m 

—  x  10' 

— 

q2 

4tt 
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Physical 

Quantity 


Dimensions 


Permittivity 

Polarization 

Potential 

Power 

j 

Power 

density 

Pressure 

, 

Reluctance 

Resistance 

Resistivity 

Thermal  con¬ 
ductivity 

Time 

Vector 

potential 

Velocity 

Viscosity 

Vorticity 
Wo  r  k 


Svm- 

bol 

SI 

Gaussian 

j 

e 

t~q~ 

1 

ml 3 

P 

<7 

1/2 
m  ' 

r- 

pr-t 

V,4> 

ml 2 

mt/2  f  1/2 

t2q 

t 

1 

P 

ml 2 

ml2 

t3 

t3 

m 

m 

It 3 

It3 

P,P 

m 

TJX 

It3 

it 2 

71 

O 

r 

l 

ml 2 

T 

R 

ml2 

t 

iq 2 

! i 

ml 3 

t 

1,  P 

ro 

k  k 

ml 

ml 

t3 

t3 

t 

t 

t 

A 

ml 

ml'2l1''2 

tq 

t 

l 

l 

V 

— 

-* 

t 

t 

m 

m 

V- M 

— 

— 

It 

It 

I 

1 

C 

— 

— 

t 

t 

w 

ml 2 

ml2 

t2 

t2 

SI 

Units 

Conversion 

Factor 

Gaussian 

U  nits 

farad/m 

36rr  x  109 

coulomb/m” 

3  x  105 

statcouiomb 
/  cm2 

volt 

1 

-  x  10  2 

3 

statvolt 

watt 

107 

erg/sec 

watt/m3 

10 

erg/ cm2-sec 

pascal 

10 

dyne/cm2 

ampere-turn 

/weber 

4x  x  10-9 

cm-1 

ohm 

1  _ i  j 

-  x  10  11 

9 

sec/cm 

ohm-m 

1  -9 
-  x  10  9 

9 

sec 

watt/ m- 
deg  (K) 

105 

l 

erg/cm-sec- 
deg  (K) 

second  (s) 

1 

second  (sec) 

weber/m 

106 

gauss-cm 

m/s 

102 

cm/sec 

l 

kg/m-s 

10 

poise 

1 

s-1 

1 

sec- 1 

joule 

107 

erg 
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INTERNATIONAL  SYSTEM  (SI)  NOMENCLATURE5 


Physical 

Name 

Physical 

Name 

Symbol 

Quantity 

of  Unit 

Quantity 

of  Unit 

for  Unit 

•length 

meter 

m 

electric 

potential 

volt 

V 

*mass 

kilogram 

kg 

electric 

ohm 

a 

•time 

second 

s 

resistance 

•current 

ampere 

A 

electric 

conductance 

siemens 

s 

•temperature 

kelvin 

K 

electric 

farad 

F 

•amount  of 

mole 

mol 

capacitance 

substance 

magnetic  flux 

weber 

Wb 

•luminous 

candela 

cd 

intensity 

magnetic 

inductance 

henry 

H 

jplane  angle 

radian 

rad 

magnetic 

tesla 

T 

fsolid  angle 

steradian 

sr 

intensity 

frequency 

hertz 

Hz 

luminous  flux 

lumen 

111! 

energy 

joule 

J 

illuminance 

lux 

lx 

force 

newton 

N 

activity  (of  a 
radioactive 

becquerel 

Bq 

pressure 

pascal 

Pa 

source) 

power 

watt 

W 

absorbed  dose 
(of  ionizing 

gray 

Gy 

electric  charge 

coulomb 

C 

radiation) 

•SI  base  unit  fSupplementary  Unit 


METRIC  PREFIXES 


Multiple 

Prefix 

Symbol 

Multiple 

Prefix 

Symbol 

10~l 

deci 

10 

deca 

da 

W 

1 

O 

H 

centi 

102 

hecto 

h 

icr3 

milli 

103 

kilo 

k 

icr6 

micro 

106 

mega 

M 

10-9 

nano 

109 

giga 

G 

io-12 

pico 

P 

10 12 

tera 

T 

10-15 

femto 

1015 

peta 

P 

Iq-13 

atto 

mmmt 

1018 

exa 

E 
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PHYSICAL  CONSTANTS  (SI)7 


Physical  Quantity  j 

Symbol 

Value 

Units 

Boltzmann  constant 

k 

1.3S07  x  10"23 

JIC1 

Elementary  charge 

e 

1.6022  x  10~19 

c 

Electron  mass 

TTt  £ 

9.1095  x  10~31 

kg 

Proton  mass 

T7X  p 

1.6726  x  10_2T 

kg 

Gravitational  constant 

G 

6.6720  x  10-u 

m3s-2kg“ 

Planck  constant 

h 

6.6262  x  10~34 

Js 

y. 

cu 

II 

1.0546  x  10-34 

Js 

Speed  of  light  in  vacuum 

C 

2.9979  x  108 

ms-1 

Permittivity  of 

fo 

8.S542  x  10~12 

F  m  ~ 1 

free  space 

Permeability  of 

/io 

4?r  x  10-7 

H  m-1 

free  space 

Protou/electron  mass 

mp/me 

1.8362  x  103 

ratio 

Electron  charge/mass 

e/me 

1.7588  x  1011 

Ckg"1 

ratio 

me 4 

Rydberg  constant 

1.0974  x  10' 

m  1 

Bohr  radius 

a  o  =  e0h2  fir  me" 

5.2918  x  10_u 

m 

Atomic  cross  section 

irao2 

8.7974  x  10“21 

•y 

m* 

Classical  electron  radius 

re  =  e2  /  47r  Come2 

2.8179  x  10-15 

m 

Thomson  cross  section 

(8t /3)r  e2 

6.6524  x  10-29 

*» 

in* 

Compton  wavelength  of 

h/mec 

2.4263  x  10-12 

m 

electron 

h/me  c 

3.8616  x  10-13 

m 

Fine-structure  constant 

a  =  e2/2fo he 

7.2974  x  10-3 

a"1 

137.04 

First  radiatiou  constant 

ci  =  2t/ic* 

3.741S  x  10-2 

Wm2 

Second  radiatiou 

C2  =  hc/k 

1.43S8  x  10"2 

mK 

constant 

Stefan- Boltzmann 

a 

5.6703  x  10“8 

W  m  '  2  K  “ 

constant 

t 
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Physical  Quantity 

Symbol 

Value 

Wavelength  associated 

Ao  =  hc/e 

1.2399  x  lO"6 

with  1  eV 

Frequency  associated 

i/0  =  e/h 

2.4180  x  1014 

with  1  eV 

Wave  number  associated 

© 

II 

n> 

o 

8.0655  x  10s 

with  1  eV 

Energy  associated  with 

hi/  o 

1.6022  x  10“19 

1  eV 

Energy  associated  with 

he 

1.9865  x  10"25 

1  m-1 

Energy  associated  with 

me3 /&€q7  h2 

13.606 

1  Rydberg 

Energy  associated  with 

k/e 

8.6173  x  10-5 

1  Kelvin 

Temperature  associated 

e/k 

1.1605  x  104 

with  1  eV 

Avogadro  number 

na 

6.0220  x  1023 

Faraday  constant 

F  =  NAe 

9.6485  x  104 

Gas  constant 

R  -  NAk 

8.3144 

Loschmidt’s  number 

n0 

2.6868  x  1025 

(no.  density  at  STP) 

Atomic  mass  unit 

mu 

1.6606  x  10"27 

Standard  temperature 

To 

273.16 

Atmospheric  pressure 

p  o  —  no  kTo 

1.0133  x  105 

Pressure  of  1  mm  Hg 

1.3332  x  102 

(1  torr) 

Molar  volume  at  STP 

Vo  =  RTq  /  po 

2.2415  x  10“2 

Molar  weight  of  air 

Mair 

2.8971  x  10"2 

calorie  (cal) 

4.1868 

Gravitational 

9 

9.8067 

acceleration 

Units 


mol  1 
C  mol-1 
J  K-1mol-1 
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PHYSICAL  CONSTANTS  (cgs)7 


Physical  Quantity 

Symbol 

Value 

Units 

Boltzmann  constant 

k 

1.3807  x  10-16 

erg/ deg  (K) 

Elementary  charge 

e 

4.8032  x  10“10 

statcoulomb 

Electron  mass 

TOe 

9.1095  x  10'2S 

(statcoul) 

g 

Proton  mass 

mp 

1.6726  x  10~24 

g 

Gravitational  constant 

G 

6.6720  x  10-8 

dyue-cm2/g2 

Planck  constant 

h 

6.6262  x  10~27 

erg-sec 

h  =  h/2ir 

1.0546  x  10~27 

erg-sec 

Speed  of  light  in  vacuum 

c 

2.9979  x  1010 

cm/sec 

Proton/electron  mass 

mp/me 

1.8362  x  103 

ratio 

Electron  charge/mass 

e/me 

5.2728  x  101T 

statcoul/g 

ratio 

Rydberg  constant 

n  2  ff2me4 

Ro°  ~  ch* 

1.0974  x  10s 

cm 

Bohr  radius 

ao  =  h2 /me2 

5.2918  x  10“9 

cm 

Atomic  cross  section 

irao 2 

8.7974  x  10~1T 

cm’ 

Classical  electron  radius 

re  =  e2 /me2 

2.8179  x  10~13 

cm 

Thomson  cross  section 

(8*/3)re2 

6.6524  X  10~23 

cm2 

Compton  wavelength  of 

h/m^c 

2.4263  X  10~10 

cm 

electron 

h/m^c 

3.8616  X  10~11 

cm 

Fine-structure  constant 

ot  —  e2 /he 

7.2974  X  10”3 

First  radiation  constant 

a~l 

d  =  2ir  he2 

137.04 

3.7418  x  10“5 

erg-cm2/sec 

Second  radiation 

C2  =  hc/k 

1.4388 

cm-deg  (K) 

constant 

Stefan- Boltzmann 

<T 

5.6703  x  10-5 

erg/cm2- 

constant 

Wavelength  associated 

Aq 

1.2399  x  10-4 

sec-deg4 

cm 

with  1  eV 
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Physical  Quantity 


Symbol 


Value 


Units 


Frequency  associated 

v0 

2.4180  x  1014 

with  1  eV 

Wave  number  associated 

fco 

8.0655  x  103 

with  1  eV 

Energy  associated  with 

1.6022  x  10"12 

1  eV 

Energy  associated  with 

1.9865  x  10'16 

1  cm"1 

Energy  associated  with 

13.606 

1  Rydberg 

Energy  associated  with 

8.6173  x  10"s 

1  deg  Kelvin 

Temperature  associated 

1.1605  x  104 

with  1  eV 

Avogadro  number 

Na 

6.0220  x  1023 

Faraday  constant 

F  =  N  a  e 

2.8925  x  1014 

Gas  constant 

R  =  NAk 

8.3144  x  107 

Loschmidt’s  number 

n0 

2.6868  x  1019 

(no.  density  at  STP) 

Atomic  mass  unit 

mu 

1.6606  x  10'24 

Standard  temperature 

To 

273.16 

Atmospheric  pressure 

po  =  nokTo 

1.0133  x  106 

Pressure  of  1  mm  Hg 

1.3332  x  103 

(1  torr) 

Molar  volume  at  STP 

Vo  =  RTo/po 

2.2415  x  104 

Molar  weight  of  air 

Mair 

28.971 

calorie  (cal) 

4.1868  x  107 

Gravitational 

9 

980.67 

acceleration 

Hz 

cm-1 

erg 

erg 

eV 

eV 

deg  (K) 
mol"1 

statcoul/mol 

erg/deg-mol 

cm-3 

g 

deg  (K) 
dyne/cm2 
dyne/ cm2 

cm3 

g 

erg 

cm/sec2 
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FORMULA  CONVERSION8 

Here  a  =  102  cm.  m~ 1 ,  /3  =  10'  erg  3  “ 1 ,  eo  —  8.8542xl0-l2F  m-v , 
Hq  =  4rr  x  1 0  ~  ‘  Hm"1,  c  =  (e0Mo)"1/2  =  2.9979  x  108ms"1,  and 
h  =  1.0546  x  10-34  J  s.  To  derive  a  dimensionally  correct  SI  formula 
from  one  expressed  in  Gaussian  units,  substitute  for  each  quantity  ac¬ 
cording  to  Q  —  kQ,  where  k  is  the  coefficient  in  the  second  column 
of  the  table  corresponding  to  Q  (overbars  denote  variables  expressed  in 
Gaussian  units).  Thus,  the  formula  do  =  h~ /fixe2  for  the  Bohr  radius 
becomes  aa0  =  (h0)2 /[(m0 /a2)(e2  oc0 /4ireo)],  or  <*o  =  eoh2 /irme2 .  To 
go  from  SI  to  natural  units  in  which  h  =  c  =  1  (distinguished  by  a  cir¬ 
cumflex),  use  Q  —  k~lQ,  where  k  is  the  coefficient  corresponding  to  Q 
in  the  third  column.  Thus  do  =  Aireoh2 /[(mh/c)(e2  eohc)]  =  47r/me2. 
(In  transforming  from  SI  units,  do  not  substitute  for  eo,  Mo,  or  c.) 


Physical  Quantity 

Gaussian  Units  to  SI 

Natural  Units  to  SI 

Capacitance 

a/4ireo 

eo-1 

Charge 

(a0/ 4ire0)1/2 

(e0dc)'1/2 

Charge  density 

(0/4vctseo)l/2 

(eo^c)“1/2 

Current 

( a0/4ireo)l,t 

(m/Ac)1'1 

Current  density 

(/3/4tt  a3«o)l/2 

(mo/M1/2 

Electric  field 

(4T/?€0/a3)1/2 

(«o /hc)1/2 

Electric  potential 

(4tt0€o/ a)l/2 

(e0/*c)l/2 

Electric  conductivity 

(4Te0)"1 

eo"1 

Energy 

P 

(Ac)-1 

Energy  density 

P/a2 

(Ac)"1 

Force 

P/a 

(Ac)'1 

Frequency 

1 

c  1 

Inductance 

4-xf.o/a 

Mo"1 

Length 

a 

1 

Magnetic  induction 

(4*P  /  a2  no)1'2 

(/i0dc)"1/2 

Magnetic  intensity 

(4*^0 /a2)1'2 

(m/Ac)1/2 

Mass 

P/a 2 

c/h 

Momentum 

P/a 

h~l 

Power 

0 

(he2)-1 

Pressure 

P/a 3 

(he)-1 

Resistance 

47T€o /a 

(«o/Mo)l/2 

Time 

1 

c 

Velocity 

a 

c~ 1 
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MAXWELL’S  EQUATIONS 


Name  or  Description 

SI 

Gaussian 

Faraday's  law 

5B 

V  x  E  =  -  — 

dt 

_  _  1  dB 

V  x  E  = - 

c  dt 

Ampere’s  law 

3D 

V  x  H  =  — -  +  J 

dt 

1  3D  4t 

V  x  H  =  -—  +  — J 
c  dt  c 

Poisson  equation 

V  •  D  =  p 

V  •  D  =  47rp 

[Absence  of  magnetic 
monopoles] 

V  •  B  =  0 

V  •  B  =  0 

Lorentz  force  on 

q  (E  +  v  x  B) 

q  (  E  +  -v  x  B  ] 

charge  q 

\  c  / 

Constitutive 

D  =  eE 

D  =  fE 

relations 

B  =  /zH 

B  =  /zH 

In  a  plasma,  /z  ss  p.0  =  4ir  X  lO""'  Hm  1  (Gaussian  units:  p.  «  1). 
The  permittivity  satisfies  €  «  €q  —  8.8542  x  10“ 12  F  m_1  (Gaussian: 
<  W  1)  provided  that  all  charge  is  regarded  as  free.  Using  the  drift 
approximation  vx  =  E  X  B/B 2  to  calculate  polarization  charge  density 
gives  rise  to  a  dielectric  constant  K  =  f/eo  =  1  +  36*  x  10 0 p/ B2  (SI) 
=  1  +  4tt pc2 / B2  (Gaussian),  where  p  is  the  mass  density. 

The  electromagnetic  energy  in  volume  V  is  given  by 


W=^j  dV(H  •  B  +  E  •  D)  (SI) 

=  —  /  dV{  H  ■  B  +  E  •  D)  (Gaussian). 


Poynting's  theorem  is 


d\V 

dt 


+  /  N  •  dS  =  -  /  dV J  •  E, 


where  S  is  the  closed  surface  bounding  V  and  the  Poynting  vector  (energy 
flux  across  5)  is  given  by  N  =  E  x  H  (SI)  or  N  =  cE  x  H/47T  (Gaussian). 


19 


ELECTRICITY  AND  MAGNETISM 


In  the  following,  t  =  dielectric  permittivity.  ;r  =  permeability  of 
conductor,  ft '  —  permeability  of  surrounding  medium,  a  =  conductivity. 
/  =  / 2tr  —  radiation  frequency,  —  /r/^o  and  n.e  =  f/fo-  Whore  sub¬ 

scripts  arc  used.  ‘1’  denotes  a  conducting  medium  and  ‘2'  a  propagating 
(lossless  dielectric)  medium.  All  units  are  SI  unless  otherwise  specified. 


Permittivity  of  free  space 

Permeability  of  free  space 

Resistance  of  free  space 

Capacity  of  parallel  plates  of  area 
A .  separated  by  distance  d 

Caj  acity  of  concentric  cylinders 
of  length  1,  radii  a,  b 

Capacity  of  concentric  spheres  of 
radii  a .  b 

Self-inductance  of  wire  of  length 
l ,  carrying  uniform  current 

Mutual  inductance  of  parallel  wires 
of  length  l ,  radius  a,  separated 
by  distance  d 

Inductance  of  circular  loop  of  radius 
6,  made  of  wire  of  radius  a, 
carrying  uniform  current 

Relaxation  time  in  a  lossy  medium 

Skin  depth  in  a  lossy  medium 

Wave  impedance  in  a  lossy  medium 

Transmission  coefficient  at 
conducting  surface 
(good  only  for  T  <g.  1)° 

Field  at  distance  r  from  straight  wire 
carrying  curreut  I  (amperes) 

Field  at  distance  z  along  axis  from 
circular  loop  of  radius  a 
carrying  current  I 


( o  =  S.S542  x  10-12  Fm'1 

Ho  =  4t  x  10~‘  Hm"1 
=  1.25G6  x  10“c  H  m_l 

R0  =  (mo  /fo)X/2  =  37G.73Q 
C  =  eA/d 

C  =  2  trel\u(b/ a) 

C  =  4: reab/(b  —  a) 

L  ~  fd 

L  =  (h'I/Att)  (1  +  4  ln(d/a)] 

L  =  &{/i'[lu(8&/a)  -2]  +  /t/4} 

r  —  e/<7 

6  =  (2 /v(U7)l/2  =  (Trftur)"'1'2 

Z  =  [;t/(€  -F  icr/w))1/2 
T  =  4.22  X  10“4(/KmlKe2/<7)l/2 

Ds  =  ^I/2-xr  tesla 

=  0.21/ r  gauss  (r  in  cm) 

B  z  =  Ma2r/[2(a3  +  z2)3/2] 
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ELECTROMAGNETIC  FREQUENCY/ 
WAVELENGTH  BANDS10 


Designation 

Frequency  Range 

Wave  [eng 

th  Range 

Lower 

Upper 

Lower 

Upper 

ULF* 

10  Hz 

3  Mm 

ELF* 

10  Hz 

3kHz 

100  km 

3  Mm 

VLF 

3  kHz 

30  kHz 

10  km 

100  km 

LF 

30  kHz 

300kHz 

1  km 

10  km 

MF 

300  kHz 

3  MHz 

100  m 

1  km 

HF 

3  MHz 

30  MHz 

10  m 

100  m 

VHF 

30  MHz 

300  MHz 

1  m 

10  m 

UHF 

300  MHz 

3  GHz 

10  cm 

1  m 

SHF  t 

3  GHz 

30  GHz 

1  cm 

10  cm 

S 

2.6 

3.95 

7.6 

11.5 

G 

3.95 

5.85 

5.1 

7.6 

J 

5.3 

8.2 

3.7 

5.7 

H 

7.05 

10.0 

3.0 

4.25 

X 

8.2 

12.4 

2.4 

3.7 

M 

10.0 

15.0 

2.0 

3.0 

P 

12.4 

18.0 

1.67 

2.4 

K 

18.0 

26.5 

1.1 

1.67 

R 

26.5 

40.0 

0.75 

1.1 

EHF 

30  GHz 

300  GHz 

1  mm 

1  cm 

Submillimeter 

300  GHz 

3  THz 

100  fim 

1  mm 

Infrared 

3  THz 

430  THz 

700  nm 

100  /zm 

Visible 

430  THz 

750  THz 

400  nm 

700  nm 

Ultraviolet 

750  THz 

30  PHz 

10  uni 

400  nm 

X  Ray 

30  PHz 

3  EHz 

100  pm 

10  nm 

Gamma  Ray 

3  EHz 

100  pm 

Note:  In  spectroscopy  the  angstrom  (A)  is  sometimes  used  (1  A  =  10~8  cm 
=  0.1  nm). 

*The  boundary  between  ULF  and  ELF  is  variously  defined. 

|The  SHF  (microwave)  band  is  further  subdivided  approximately  as  shown. 
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AC  CIRCUITS 


For  a  resistance  R.  inductance  L.  and  capacitance  C  in  series  with 
a  voltage  source  V  =  V0  exp(iwt)  (here  t  =  \f  —  1),  the  current  is  given 
by  /  =  dq/dt ,  where  q  satisfies 


dt 


=  V. 


Solutions  are  q(t )  =  q,  +  qt .  /(t)  =  I,  +  It.  where  the  steady  state  is 
/,  =  iu/q,  =  V/Z  in  terms  of  the  impedance  Z  =  R  +  i(u>L  —  1/uC)  and 
It  —  dqt/dt.  For  initial  conditions  <?(0)  =  qo  =  5o  +  qj,  I( 0)  S  7o,  the 
transients  can  be  of  three  types,  depending  on  A  =  il"  —  4 L/C: 

(a)  Overdamped,  A  >  0 


A 


7p  4-  7+  q  o 
7+  ~  7- 


exp(-7-t) 


7p  -f  7-  go 
7+  -  7- 


exp(-7+t), 


7+(7a+7-<?o)  ,  7-(^o+7+9o) 

- exp(— 7+t) - exp(  — 7_t), 

7+  “  7-  7+  ~  7- 


where  7i  =  (R  ±  Al^2)/2L; 

(b)  Critically  damped,  A  =  0 

qt  =  (qo  +  (fo  +  7Rqo)t]  exp(-7Rt), 
ft  =  [7o  -  (7o  +  7R7o)7Rt]  exp(-7«t), 


where  7 r  =  R/2L-, 

(c)  Underdamped,  A  <  0 


qe 

It 


7R7o  +  lo  .  ,  .  _  .1  ,  . , 

- sin  t  +  qo  cos  u\ t  exp(— 7Rt), 

w; 

r  „  (wt2  +  7RZ)-7o  +  7r/o  .  , 

/ o  coswjt  —  - sin(wit) 

Wi 


cxp(  -7r£), 


where  =  w0(l  -  R2C/4L)l/2  and  w0  =  (IC)_l/2  is  the  resonant 
frequency.  At  u  =  w0,  Z  =  R.  The  quality  of  the  circuit  is  Q  =  uqL/R. 
Instability  results  when  L,  R ,  C  are  not  all  of  the  same  sign. 
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DIMENSIONLESS  NUMBERS  OF  FLUID  MECHANICS12 


Namei  s ) 

Symbol 

Definition 

Significance 

Alfven, 

Karman 

Al.  Ka 

vA/v 

*(Magnetic  force/ 
inertial  force)1/2 

Bond 

Bd 

{p  -p)L2g/X 

Gravitational  force/ 
surface  tension 

Boussinesq 

B 

V7(2  gR)1/2 

(Inertial  force/ 

gravitational  force)1/2 

Brinkman 

Br 

nV2/kAT 

Viscous  heat/conducted  heat 

Capillary 

Cp 

pV/Yl 

Viscous  force/surface  tension 

Carnot 

Ca 

(T2  -to/t2 

Theoretical  Carnot  cycle 
efficiency 

Cauchy, 

Hooke 

Cy,  Hk 

pv2/ r  =  m2 

Inertial  force/ 

compressibility  force' 

Clausius 

Cl 

LV3  p/kAT 

Kinetic  energy  flow  rate/heat 
conduction  rate 

Cowling 

C 

(VA/V)2  =  Al2 

Magnetic  force/inertial  force 

Crispation 

Cr 

P-k/HL 

Effect  of  diffusion/effect  of 
surface  tension 

Dean 

D 

D3,2V/v{2t)x/2 

Transverse  flow  due  to 

curvature/lougitudiual  flow 

[Drag 

coefficient] 

cD 

(p'  -  p)Lg / 
p'v 2 

Drag  force/iuertial  force 

Eckert 

E 

V2  /  CpAT 

Kinetic  energy/change  in 
thermal  energy 

Ekman 

Ek 

(u/2CIL2)1/2  = 
(Ro/Re)1/2 

(Viscous  force/Coriolis  force)1/2 

Euler 

Eu 

Ap/pV 2 

Pressure  drop  due  to  friction/ 
dynamic  pressure 

Fronde 

Fr 

v/(aL)1/2 

V/NL 

f(Inertial  force/gravitatioual  or 
buoyancy  force)1/2 

Gay-Lussac 

Gti 

1  /PAT 

Inverse  of  relative  change  in 
volume  during  heating 

Grashof 

Gr 

gL3pAT/u2 

Buoyancy  force/viscous  force 

[Hall 

coefficient] 

Ch 

Gyrofreqttency/ 
collision  frequency 

*(t)  Also  defined  as  the  inverse  (square)  of  the  quantity  shown. 
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Name(s)  j 

Symbol 

Definition  ) 

Significance 

Hartmann 

H 

BL/ifi  V)1/2  = 
(Rm  Re  C) 1/2 

Magnetic  force/ 
dissipative  force 

Kuudsen 

Kn 

X/L 

Hydrodynamic  time/ 
collision  time 

Lorentz 

Lo 

V/c 

Magnitude  of  relativistic  effects 

Lundquist 

Lu 

PqLVa  / t)  - 
A1  Rm 

J  x  B  force/resistive  magnetic 
diffusion  force 

Mach 

M 

V/Cs 

Magnitude  of  compressibility 
effects 

Magnetic 

M  ach 

Mm 

V/VA  =  Al"1 

(Inertial  force/magnetic  force)1/2 

Magnetic 

Reynolds 

Rm 

Mo  LV/V 

Flow  velocity/magnetic  diffusion 
velocity 

Newton 

Nt 

F/pL2V 2 

Imposed  force/inertial  force 

Nusselt 

N 

aL/k 

Total  heat  transfer/thermal 
conduction 

Peclet 

Pe 

LV/k 

Heat  convection/heat  conduction 

Poisseuille 

Po 

D2±p/pLV 

Pressure  force/viscous  force 

Prandtl, 

Schmidt 

Pr,  Sc 

v/k 

Momentum  diffusion/ 
heat  diffusion 

Rayleigh 

Ra 

gH3pAT/uK 

Buoyancy  force/diffusion  force 

Reynolds 

Re 

LV/u 

Inertial  force/viscous  force 

Richardson 

Ri 

( NH/AV )2 

Buoyancy  effects/ 
vertical  shear  effects 

Rossby 

Ro 

V/2QL  sin  A 

Inertial  force/Coriolis  force 

Stanton 

St 

a/pcpV 

Thermal  conduction  loss/ 
heat  capacity 

Stefan 

Sf 

a  LT3  /  k 

Radiated  heat/couducted  heat 

Stokes 

S 

u/L2f 

Viscous  damping  rate/ 

■vibration  frequency 

Strouhal 

Sr 

fL/V 

Vibration  speed/flow  velocity 

Taylor 

Ta 

(2QL2/u)2 

Rl/2{AR)3/2 

■{Sl/v) 

Centrifugal  force/viscous  force 
(Centrifugal  force/ 
viscous  force)1/2 

Thring. 

Boltzmann 

Th.  Bo 

pcpV/€aT3 

Convective  heat  transport/ 
radiative  heat  transport 

We  b  e  r 

W 

pLV2 /E 

Inertial  force/surface  tension 
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Nomenclature: 


B 

C , ,  c 

cp 

D  —  2R 

F 

f 

9 

H,  L 

k  —  pCpK 

N  =  (g/H)1/2 

R 

r 

tl 

T 

V 

VA  =  B/(Mop)1/2 
a 

(3 

r 


e 

V 

K 

A 

A 

p  =  pu 

Mo 

v 

P 

t 

P 

E 

a 

a 


Magnetic  induction 
Speeds  of  sound,  light 

Specific  heat  at  constant  pressure  (units  m"  s~2  K-1) 

Pipe  diameter 

Imposed  force 

Vibration  frequency 

Gravitational  acceleration 

Vertical,  horizontal  length  scales 

Thermal  conductivity  (units  kgm-1s~2) 

Brunt-Vaisala  frequency 

Radius  of  pipe  or  channel 

Radius  of  curvature  of  pipe  or  channel 

Larmor  radius 

Temperature 

Characteristic  flow  velocity 
Alfven  speed 

dT 

Newton’s-law  heat  coefficient,  k——  —  a£T 

dx 

Volumetric  expansion  coefficient.  dV JV  =  (3dT 

Bulk  modulus  (units  kgm-1s~2) 

Imposed  difference  in  two  radii,  velocities, 
pressures,  or  temperatures 

Surface  emissivity 

Electrical  resistivity 

Thermal  diffusivity  (units  m2  s_l) 

Latitude  of  point  on  earth’s  surface 
Collisioual  mean  free  path 
Bulk  viscosity 
Permeability  of  free  space 
Kinematic  viscosity  (units  m3  s_l) 

Mass  density  of  fluid  medium 

Mass  density  of  bubble,  droplet,  or  moving  object 
Surface  tension  (units  kgs-2) 

Stefan-Boltzmann  constant 
Solid-body  rotational  angular  velocity 
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SHOCKS 


At  a  shock  front  propagating  in  a  magnetized  fluid  at  an  angle  9 
with  respect  to  the  magnetic  induction  B.  the  jump  conditions  are  ‘3*4 

(1)  pU  =  pU  =  q: 

(2)  pU2  +  p  +  B^/2p  =  pU2  +  p  -f  Brf  2/x: 

(3)  pUV  -  S„  Bj_/ p.  =  pUV  -  B^/p; 

(4)  B,|  =  B„; 

(5)  U Bj_  -  VB,|  =  D'Bx  -  V’Bu; 

(6)  i({72  +  V2)  +  w  +  (U B ±  -  VB„  B±)/ppU 

=  §(B2  +  V2)  +  u,  +  (i/B,2  -  V  Bn  B  ±)  /  ppU . 

Here  U  and  V  are  components  of  the  fluid  velocity  normal  and  tangential 
to  the  front  in  the  shock  frame;  p  —  l/v  is  the  mass  density;  p  is  the 
pressure;  B±  —  B  sin  9,  Bj)  =  B  cos  B;  p  is  the  magnetic  permeability 
(p  =  4tt  in  cgs  units);  and  the  specific  enthalpy  is  w  =  e  +  pv,  where 
the  specific  internal  energy  e  satisfies  de  =  Tds  —  pdv  in  terms  of  the 
temperature  T  and  the  specific  entropy  a.  Quantities  in  the  region  behind 
(downstream  from)  the  front  are  distinguished  by  a  bar.  If  B  =  0,  then15 

(T)  U-U  =  [(p-p)(v-v)}1'2- 

(8)  (p  -  p){v  -  v)~x  -■  q2\ 

(9)  w  -  w  =  \  (p  -  p)(v  +  v); 

(10)  e  -  e  =  \{p  +  p)(v  -  v). 

In  what  follows  we  assume  that  the  fluid  is  a  perfect  gas  with  adiabatic 
index  7  =  1  4-  2/n,  where  n  is  the  number  of  degrees  of  freedom.  Then 
p  =  pRT/m,  where  R  is  the  universal  gas  constant  and  m  is  the  molar 
weight;  the  sound  speed  is  given  by  C,'  =  (c>p/<9p)j  —  1VV\  an<^  w  — 
7 e  =  ~!Pv/(~(  -f  1).  For  a  general  oblique  shock  in  a  perfect  gas  the 
quantity  X  =  r  ~ 1  ( V / VA  ) 2  satisfies14 

(11)  (X-i3/*){X -cos2  9)2  =  X  sin2  9  {[1  +  (r  -  l)/2a]  X  -  cos'  fl}, 
where  r  =  p/p.  a  =  ^  [  7  -I-  1  —  (7  —  l)rj,  and  p  =  C,2  jV.\  2  ~  4t7  p/ B2. 
The  density  ratio  is  bounded  by 

(12)  1  <  r  <  (7  +  1 ) / ( 7  -  1). 

If  the  shock  is  normal  to  B  (he.,  if  9  —  ir/2),  then 

(13)  U2  =  (r  /  a)  {C,2  +  V^2  (1  4-  (1  -  7/2)(r  -  1)1  } ; 
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7' 


(14)  U/U  =  B/B  - 

(15)  V  =  V; 

(16)  p  =  p  -i-  (1  -  r~l)pU2  -f  (1  -  r2)B2  /2p. 

If  8  =  0,  there  are  two  possibilities:  switch-on  shocks,  which  require 
P  <  1  and  for  which 

(17)  U2  =  rVA2; 

(18)  U  =  Va2 /U\ 

(19)  fl12  =  2B112(r-l)(Q-/3); 

(20)  V  =  UB±/Bl 

(21)  p  =  p  +  pU2(l  -  a  +  P)(l  -  r-1), 

and  acoustic  (hydrodynamic)  shocks,  for  which 

(22)  U2  =  (r/a)C.2; 

(23)  U  =  U/r; 

(24)  V  =  B±  =  0; 

(25)  p-p  +  pU2{  1  -  r-1). 

For  acoustic  shocks  the  specific  volume  and  pressure  are  related  by 

(26)  v/v  =  [(7  +  1  )p  +  (7  -  1)?]  /  [(7  -  1  )P  +  (7  +  1)?]- 
In  terms  of  the  upstream  Mach  number  M  —  U/C, , 

(27)  p/p  =  v/v  =  U/U  =  (7  +  l)M2/[(7  -  1)M2  +  2]; 

(28)  p/p  =  (2 7M2  -  7  +  l)/(7  +  1); 

(29)  T/T  =  [(7  -  1)M2  +  2](27M2  -  7  +  l)/(7  +  1)2M2; 

(30)  M2  =  ((7  -  1)M 2  +  2]/[27M2  -7  +  1]- 
The  entropy  change  across  the  shock  is 

(31)  As  =  s  —  3  =  cu  ln[(p/p)(p/p)7], 

where  cv  =  R/(~f  —  l)m  is  the  specific  heat  at  constant  volume;  here  R 
is  the  gas  constant.  In  the  weak-shock  limit  (M  — *  1), 


(32)  As  — ►  cv 


27(7  ~  1) 
3(7  +  1) 


-l)3  * 


167  R 
3(7  +  l)17* 


(M  -  1)' 


The  radius  at  time  t  of  a  strong  spherical  blast  wave  resulting  from  the 
explosive  release  of  energy  E  in  a  medium  with  uniform  density  p  is 


(33)  RS  =  C0(£?i2/p)1/s, 


where  Co  is  a  constant  depending  on  7.  For  7  =  7/5,  Co  =  1.033. 
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FUNDAMENTAL  PLASMA  PARAMETERS 


All  quantities  are  in  Gaussian  cgs  units  except  temperature  ( T .  Te. 
T, )  expressed  in  eV  and  ion  mass  [m,j  expressed  in  units  of  the  proton 
mass,  p.  =  m,/mp;  Z  is  charge  state:  k  is  Boltzmann's  constant:  K  is 
wavelength:  7  is  the  adiabatic  index:  In  A  is  the  Coulomb  logarithm. 


Frequencies 

electron  gyrofrequency 

ion  gyrofrequency 
electron  plasma  frequency 


ion  plasma  frequency 


electron  trapping  rate 
ion  trapping  rate 


electron  collision  rate 
ion  collision  rate 

Lengths 

electron  deBroglie  length 

classical  distance  of 
minimum  approach 

electron  gyroradius 
ion  gyroradius 

plasma  skin  depth 
Debye  length 


}ct  =  w„/2t  =  2. SO  x  106S  Hz 

wce  —  eB/mtc  =  1.76  x  10'  B  rad/sec 

fc,  =  oc,/2x  =  1.52  x  103Zm-1£  Hz 

u/cl  =  eB/m,c  =  0.58  x  10 3  Zfi~l  B  rad /sec 

fPe  =  wp,/2i  =  8.98  x  103nt1/2  Hz 

w,«  =  (4xn«e2/me)1/2 

=  5.64  X  104ne1^2  rad/sec 

/pi  =  Wpi  /2tr 

=  2.10  x  102Z/i-l/2n,1/2  Hz 
u»?,-  =  (4tti,  Z2e2/m,  )1/2 

=  1.32  X  105ZM_1/2n,1/2rad/sec 
t/Tc  =  (eKE/m,)1'2 

=  7.26  x  108/Cl/2£l/2  sec-1 
i/r,  =  (cK E/m,  )l/2 

=  1.69  x  107fCl/2£l/2M~l/2sec"1 
=  2.91  x  10~6ne  In  A Te~3/!  sec-1 
i/<  =  4.78  x  10"*n,  Z2  In  AT,~3/2  sec-1 

A  =  h/(mtkTt)1/2  =  2.76  x  10"*Tt_l/2  cm 
e2  /fcT  =  1.44  x  10'TT_l  cm 


=  «t,/  u/c  =  2.38Te 

1/2 

B  “ 1  cm 

r. 

=  UTt/wci 

=  1.02  x  102m1/2^* 

1 T. 

c/ 

upe  =  5.31  x  10sne" 

1/2 

cm 

Ad  =  {kT/Airne2)1'2 

=  7.43  x  102Tl/2n 

-1/ 

3  cm 
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Velocities 

electron  thermal  velocity  vre  =  (kTe  /me)1>/2 

=  4.19  x  10'  Te  1/<2  cm /sec 

ion  thermal  velocity  v-pi  =  {kTi /m,)1^2 

=  9.79  x  105/x-1/2rt1/2  cm/sec 
ion  sound  velocity  C4  =  (7 ZkTe/m,i)1^2 

=  9.79  x  105(7ZTe/ji)1/2  cm/sec 
Alfven  velocity  =  B /(4irnj  )1/2 

=  2.18  x  1011/i“1/2ni_1/2Bcm/sec 

Dimensionless 

(electron/proton  mass  ratio)1/2  (me/rrip)1/2  =  2.33  X  10“2  =  1/42.9 

number  of  particles  in  (4-jr/3)nXD3  =  1-72  X  10°T3/2n.~1/2 

Debye  sphere 

Alfven  velocity/speed  of  light  vA/c  =  7.28/i-1/2nt-  “l/2B 

electron  plasma/gyrofrequency  wpe/wce  =  3.21  x  10-3ne1/2 B-1 
ratio 

ion  plasma/gyrofrequency  ratio  u/p,  /o»c,-  =  0.  137/i1/2n,1/2B~l 
thermal/magnetic  energy  ratio  P  =  8irnkT/ B2  =  4.03  x  10“11nTjB-2 
magnetic/ion  rest  energy  ratio  B2 /8irnj m^c2  =  26.5£i~1ni  ~1  B2 

Miscellaneous 

Bohm  diffusion  coefficient  Db  =  ( ckT/X6eB ) 

=  6.25  x  106TB~l  cm2 /sec 

transverse  Spitzer  resistivity  =  1-15  x  10~14Z  In  AT‘  ~3/2  sec 

=  1.03  x  10-2Z  In  AT” 3/2  Q  cm 

The  auomalous  collision  rate  due  to  low-frequency  ion-sound  turbulence 
is  ^  _ 

u*  «  upeW/kT  =  5.64  x  104ne1/2VV  /kT  sec-1, 

where  W  is  the  total  energy  of  waves  with  u> / K  <  vt»- 
Magnetic  pressure  is  giveu  by 

Pm»s  =s  B2  /8ir  =  3.98  x  10®  B2  dynes/cm2  =  3.93(B/B0)2  atm, 
where  Bo  =  10  kG  =  1  T. 

Detonation  energy  of  1  kilotou  of  high  explosive  is 

WjtT  =  1012  cal  =  4.2  X'1019  erg. 
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PLASMA  DISPERSION  FUNCTION 


Definition16  (first  form  valid  only  for  Im  C  >  0): 


z(  0  =  *~1/2 


dt  exp  (~^2) 


t  -  c 


=  2; exp 


(-C2)  f  dt  exp  (-r)  ■ 
J  —  oo 


Physically,  (  =  i  +  iy  is  the  ratio  of  wave  phase  velocity  to  thermal 
velocity. 

Differential  equation: 

dZ  ,  /,  d2Z  dZ 

—  =  -2(1  +  <Z)  ,  Z(0)  =  :x1/2;  —  +  2C~  +  2Z  =  0. 

Real  argument  (y  =  0): 

Z(z)  =  exp  (  —  z2)  ^iir1/2  —  2  j  dt  exp  (t2)^  • 

Imaginary  argument  (z  =  0): 

Z(iy)  =  t*1/2  exp  (?2)  [1  -  erf(y)]  . 

Power  series  (small  argument): 

Z( 0  =  »7rW2  exp  (-C2)  -  2<  (l  -  2<2/3  +  4C“/15  -  8C6/105  +  •••). 
Asymptotic  series,  (Cl  1  (Ref.  17): 

Z( C)  =  »7r1/2(T  exp  (~C2)  — C_l  (l  +  l/2C2+3/ 4C"*  +  15 /8C6  +  •  •  •)  » 
where 

f  o  y  >  i*r 

^  ~  I  i  |y|  <  M_l 

l  2  y  <  -|zi"1 

Symmetry  properties  (the  asterisk  denotes  complex  conjugation): 

z«V  =  -[*(-01*; 

Z(C)  =  [Z(0]  *  +  2i^r1/2  exp[-(C*)2]  (y  >  0). 

Two-pole  approximations18  (good  for  £  in  upper  half  plane  except  when 

y  <  irl/2z2  exp(  — z2),  z  >  1): 

,  ,  0.50  +  0.81i  0.50-0.81i 

Z( O  as  - 7 - 7~T7 — •  a  =  0.51  -  0.81»; 


a  -  C  a*  +  C 

0.50  +  0.96i  0.50  -  0.96i 


Z'( O  ft  — —  '  - —  +  6  =  0.48-0.9K. 

(b-  <)2  (**  +  C)2 
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COLLISIONS  AND  TRANSPORT 


Temperatures  are  in  eV;  the  corresponding  value  of  Boltzmann’s  constant 
is  k  =  1.60  x  10-12  erg/eV:  masses  fis  /i'  are  in  units  of  the  proton  mass* 
ea  =  Zae  is  the  charge  of  species  a.  All  other  units  are  cgs  except  where 
noted. 

Relaxation  Rates 

Rates  are  associated  with  four  relaxation  processes  arising  from  the  in¬ 
teraction  of  test  particles  (labeled  a)  streaming  with  velocity  vQ  through  & 
background  of  field  particles  (labeled  (3): 


slowing  down 

a/(3 

—  =  v» 

transverse  diffusion 

Tt (v“  -  =  ■'I"-* 

parallel  diffusion 

|<v.  -  *.)J  - 

energy  loss 

1-y  = 

where  the  averages  are  performed  over  an  ensemble  of  test  particles  and  a 
Maxwellian  field  particle  distribution.  The  exact  formulas  may  be  written10 

v?/0  =  (l  +  rnQ/-mp)t{>{xa/0)v 

vl/0  =  2  [(1  -  l/2za/0)ifi(xa ,0)  +  Tp'(xa/I3)]  ^/0; 

»Z/(3  =  [iP(xa/0)/xa/0] 

1>Z/0  =  2  [(Tna/m3)i>(xa/0)  -  Ip'(xa/0)}  t/" /3 , 

where 


v 


a, /a 

o 


2  2 

47re0  e3  A a3n3/rn<. 


xa/0  =  m0Va7 /2kT3; 


and  \a3  =  In  Aa3  is  the  Coulomb  logarithm  (see  below).  Limiting  forms  of 
t'*,  uj_  and  t/||  are  given  in  the  following  table.  All  the  expressions  shown 
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have  uuits  cm3  sec-1.  Test  particle  energy  e  and  field  particle  temperature  T 
are  both  in  eV;  =  m.i/mp  where  mp  is  the  proton  mass:  Z  is  ion  charge 
state:  in  electron-electron  and  ion-ion  encounters,  field  particle  quantities  are 
distinguished  by  a  prime.  The  two  expressions  given  below  for  each  rate  hold 
for  very  slow  1)  and  very  fast  (x° ^ 3  1)  test  particles,  respectively. 


Electron-electron 

Ul/€'  /*€'*"'  ' 
«/«' 


Slow 


5.8  x  10-6r-3/2 


Fast 


/n./A,./  «  5.8  x  10-6T-1/2 e-1 


e / e/  /  . 

v\\  /n./A, 

Electron-ion 


«/-7  n,2TaA,i 
u‘x/i/niZ2X,i 

v\\  /n'Z 
Ion-electron 

„i'‘/neZ2\ie 

u^e/neZ2Xit 

V\\  /n*Z 


./  «  2.9  x  10 


— 6T  — 1/2  -1 


—  7.7  x  10-6<-3/2 
— *  7.7  x  10-6«-3/2 
— .  3.9  x  10-6T6-s/2 


0.23  fi3/2T~3/2 
2.5  x  10-4/i1/2T-1/2e-1 


U 

Ion-ion 


:  1.2  x  10_4/i1/2T~1/2€~l 

:  1.6  x  10-V"lT'3/2 

3.2  x  10_9M-1T-1/2«~l 

i.6  x  io~5V-lr-1/2«-1 


3.9  x  10-6e-3/2 
7.7  x  10~6e“3/2 
2.1  x  10  ~°fx-lTe~s/2 


i.7xio-V/2r3/2 

1.8  x  10-7/i-l/2e-3/2 
1.7  x  1 0”4/il/2re~5/2 


nitZ2Z'  2Xi{, 


«  6.8  x  10-8^^— -  ^1  +  T 


-3/2 


— >  9.0  x  10' 


-S  fl  ,  1  \  //2 


(;v) 


e3/2 


n;/  Z2Z'2Xiif 
i/i' 


II 


ni,Z2Z'2Xii, 


In  the  same 


s  1.4  x  io- V1/2m~1t,-1/2€-1 

— >  1.8  x  10-7/i"l/2€-3/2 

s  6.3  x  10-yi/Vlr-;/2€-1 

— *  9.0  x  10"*/*l/2/x,“lT«~5/2 
limits,  the  energy  transfer  rate  follows  from  the  identity 


v{  -  2ua  -  i/±  -  i/||, 

except  for  the  case  of  fast  electrons  or  fast  ions  scattered  by  ions,  where  the 
leading  terms  cancel.  Then  the  appropriate  forms  are 

v\,{  — ►  4.2  x  10-9n.j  Z2 Xgi 

[«-3/2M_1  -  8.9  x  10 4(/x/T)1/2«-1  exp(  — 1836/te/T)]  sec-1 
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and 


1.8  x  10 ~7nl,Z2Z'*\ii, 

f  -3/2 .1/2  ,  /  .  ,  , 


-  1.1(jj.'/T)1/~c  1  expi  -n'e/T) j 


In  general,  the  energy  transfer  rate  is  positive  for  (  >  (Qm  and  nega¬ 

tive  for  e  <  <<,*.  where  x*  =  (mj/ma)fa*/Tj  is  the  solution  of  ip‘ {x*)  = 
(m.0/m j)ip(x*).  The  ratio  (a*/Tj  is  given  for  a  number  of  specific  a,  p  in  the 
following  table: 


a/{3  i/e  e/e ,  i/i  e/p 


e/D  e/T,  e/He3  e/He4 


0.98  4.8  x  10“3  2.6  x  10~3  1.8  x  10~3  1.4  x  10~3 


When  both  species  are  near  Maxwellian,  with  T,  <  Te,  there  are  just 
two  characteristic  collision  rates.  For  Z  —  1, 

i/e  =  2.9  x  10"*nAr«"3/2sec"1; 

Vi  =  4.8  x  10"*nAT,"3/2M“l/2  sec"1. 


Temperature  Isotropization 
Isotropization  is  described  by 


IT±  =  1  dT„ 

dt  2  dt 


=  -T„), 


where,  if  A  =  Tx/T|j  -  1  >  0, 


o,  ^ x/7 e  j  n.0Aaij  _ 2 

T  ma'/2(fcr|1)3/J 


—  3  +  (A  ■+•  3) 


tan-l(A1/2) 

A1/2 


If  A  <  0,  tan“l( A1/,2)/A1//2  is  replaced  by  tanh “ 1  ( —  A)1  ^2 /( —  A)1^2 .  For 
Tx  %  T„  =  r, 

4  =  8.2  x  lO"7nAT"3/2  sec"1; 

«4  =  1.9  x  10"8nAZ2M-1/2T"3/2  sec'1. 
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Thermal  Equilibration 

If  the  components  of  a  plasma  have  different  temperatures,  but  no  rela¬ 
tive  drift,  equilibration  is  described  by 


dT\ 

dt 


=  V  *:/3(Tj  -  ra). 


j 


where 


-a/3  ,  n  —  19  (  ma  771 3  )  ^  %a  n3^aJ  -1 

i/t  —  1.8  x  10  - - - — - j_  - sec 


( maTj  +  mjTa)2/2 
For  electrons  and  ions  with  T«  as  T,  =  T,  this  implies 


£*/* / n,  =  v\'* /ne  =  3.2  x  lO'^Z* \/ cm*  sec 


-9  >r2 


'3/2 _ 3  .  -1 


Coulomb  Logarithm 

For  test  particles  of  mass  ma  and  charge  ea  =  Za  e  scattering  off  field 
particles  of  mass  mj  and  charge  ej  =  Z^e,  the  Coulomb  logarithm  is  defined 
as  A  =  In  A  =  ln(rmxx/rmio).  Here  rmin  is  the  larger  of  eae3/ma3 u2  and 
h/2mai3u,  averaged  over  both  particle  velocity  distributions,  where  mai3  = 

m,  m3 /(rna  -f  m3)  and  u  =  va  —  v3;  rm%x  =  (4t  ^  ny  e^2 /kT-, )  “ 1/2 ,  where 
the  summation  extends  over  all  species  7  for  which  u2  <  vx-,2  =  kT-j/rn y.  If 
this  inequality  cannot  be  satisfied,  or  if  either  uuca~l  <  rmax  or  u ujcj~1  < 
rm»x,  the  theory  breaks  down.  Typically  A  as  10-20.  Corrections  to  the  trans¬ 
port  coefficients  are  0(A~l);  hence  the  theory  is  good  only  to  ~  10%  and  fails 
when  A  ~  1. 

The  following  cases  are  of  particular  interest: 

(a)  Thermal  electron-electron  collisions 


A 


ee 


=  23  -  ln(ne;/2re_3/2), 
=  24  -  ln(ne1/2Te-1), 


Te  <  10  eV; 
T,  >  10  eV. 


(b)  Electron-ion  collisions 


Ae.  —  A,e  =  23  -  In  (ne1/2zr;3/j), 

=  24  -  In  (nel/2Te-1), 

=  30  -  In  (n,l/2r,  _3/2Z2/i~l), 


T,me/m,  <  Te  <  10 Z‘eV; 
T,mt/m.t  <  10 Z2  eV  <  Te 
T-e  <  T,  Zme/mt. 
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(c)  Mixed  ion-ion  collisions 


A.,,  =  A,/,  =23 -In 


+  A*') 

1 - 

t: 

Ji 

1 

■ft 

l  T,  Tx,  ) 

(d)  Counterstreaming  ions  (relative  velocity  vp  =  poc)  in  the  presence  of 
warm  electrons,  kT,  /m, ,  kT,i  /mti  <vq'<  kTt/me 


A  (1/  =  A,/,  =  35  —  In 


pil'pD2  \Tt 


Fokker-Planck  Equation 
Df°  _  df 


Dt  -  at  +V.vr  +  r.v./-- 


ar 


■)  . 

*  coll 


where  F  is  an  external  force  field.  The  general  form  of  the  collision  integral  is 

or/dt)con  =  -  vv .  with 


Ja/3  =  2tt  A 


or  /3 


Ca*t32  f  j3  ,,  2, 

-  /  a  v  (u  I  —  uu)u 

/ 

•{* 


-3 


/°(v)Vv,/V')  -  —  /P(v')Vvr(v) 

TTlnr 


r#/  '> 


(Landau  form)  where  u  =  v'  —  v  and  /  is  the  unit  dyad,  or  alternatively, 


} 


2  2  /  \ 
|/*(v)VvH(v)  -  jVv  ■  [f  (v)VvVvG(v)]  1  , 


where  the  Rosenbluth  potentials  are 


G{v)  =  /  /,3(v')u d3i/ 


tf(v)=(l  +  —  )  I .  f3{v')u~l  d\' 

mg 
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If  species  a  is  a  weak  beam  (number  and  energy  density  small  compared  with 
background)  streaming  through  a  Maxwellian  plasma,  then 


>/3  _ 


a  /  3 


V  /  —  —  V 

J  o 


a  /  3 


?«/"  +  -  ►.■'V'fv/ 


a/3 


B-G-K  Collision  Operator 

For  distribution  functions  with  no  large  gradients  in  velocity  space,  the 
Fokker-Planck  collision  terms  can  be  approximated  according  to 


£U 

Dt 


~  ft)  +  I 'ti(Ft  ~  ft); 


^  =  -..(ft  -  /.)  +  -..(ft  -/.)• 

The  respective  slowing-down  rates  given  in  the  Relaxation  Rate  section 

above  can  be  used  for  va$,  assuming  slow  ions  and  fast  electrons,  with  e  re¬ 
placed  by  Ta.  (For  vee  aud  i/;;,  one  can  equally  well  use  and  the  result 
is  insensitive  to  whether  the  slow-  or  fast-test-particle  limit  is  employed.)  The 
Maxwellians  Fa  and  Fa  are  given  by 


Fa  —  na 


exp 


f  ^a(v  ~  va)2  1  \  . 

I  L  2kTa  J  ; 


where  nQ ,  va  and  Ta  are  the  number  density,  mean  drift  velocity,  and  effective 
temperature  obtained  by  taking  moments  of  fa.  Some  latitude  in  the  definition 
of  Ta  and  va  is  possible;20  one  choice  is  Te  =  T;,  T*  =  Te,  ve  =  Vj,  v,  =  ve. 

Transport  Coefficients 

Transport  equations  for  a  multispecies  plasma: 


dana 

dt 


+  nQ  V  • 


v0  =  0; 


TtXfy  Tljj 


dt 


-  V  ■  Pa  +  Za  en„ 


1 

E  +  -  v 
c 


or 


x  B 


+  R<*; 
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+PoV  ■  va  =  -V  •  qa  -  Pa  :  Vva  +  Qa. 

2  dt 

Here  da  / dt  =  d/dt  -f  va  •  V:  pa  =  nakTa,  where  k  is  Boltzmann's  constant; 
Ra  =  ^Raj  and  Qa  =  wiiere  atl<J  Q°3  arc  respectively 

the  momentum  and  energy  gained  by  the  ath  species  through  collisions  with 
the  ,5th:  Pa  is  the  stress  tensor;  and  qQ  is  the  heat  flow. 

The  transport  coefficients  in  a  simple  two-component  plasma  (electrons 
and  singly  charged  ions)  are  tabulated  below.  Here  ||  and  _L  refer  to  the  di¬ 
rection  of  the  magnetic  field  B  =  bB;  u  =  ve  —  v;  is  the  relative  streaming 
velocity;  ne  =  n<  =  n:  j  =  — neu  is  the  current;  ujce  =  1.76  x  10'  B  sec-1  and 
u>ci  =  (me/m,)cjCe  are  the  electron  and  ion  gyrofrequencies,  respectively;  and 
the  basic  collisional  times  are  taken  to  be 


3^rn:(kTe)3/2 
4\/2tc  nAe4 


T  3/2 

=  3.44  x  10s  — - - sec, 

n  A 


where  A  is  the  Coulomb  logarithm,  and 


4v/7rnAe4  nA 


In  the  limit  of  large  fields  ( uicara  1,  a  =  t,  e)  the  transport  processes  may 

be  summarized  as  follows:21 


momentum  transfer 

Rr,= 

— R,e  =  R  =  Ru  +  Rr; 

frictional  force 

Ru  = 

ne(j|t/°'ll  +  J-l/o-jl); 

electrical 

conductivities 

*\  1  = 

ne2re 

2.0<7X  =  2.0 - 

T7X  ^ 

thermal  force 

Rr  = 

3  Tl 

-0.71nV,,(fcre) - b  x  Vx(l kTe); 

2  wcere 

ion  heating 

Qi  = 

3m,  nk 

- (Te  -  T,); 

m.i  re 

electron  heating 

Q*  = 

—  Qi  -  R  •  u; 

ion  heat  flux 

qi  = 

-xjVuCir.)  -  /civx(fcT,)  +  K\b  x  Vx(fcTi); 

ion  thermal 
conductivities 

"i  = 

nfcT.r,  ■  2nkT,  ,  5nfcr, 

3.9  t  ^  J_  2  ’  r,  ' 

m,  miUc~Tt  Zm,ujc, 

electron  heat  flux 

qe  = 

e  ,  e 

qu  +  qr; 

frictional  heat  flux 

n 

%>  3 

c r 

3nkTe  , 

0.71nkTeuu  #-f  - b  x  u  x ; 

2wceTe 
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thermal  gradient 
heat  flux 

electron  thermal 
conductivities 

stress  tensor  (both 
species) 


qr  =  -*Ji vN(fcr«)  -  v,(trt)  -  *eAb  x  vx(fcre): 

5  nkTt 


..  e  _  „  r,n^'PeTe  .  e  _  ,  „ 

^Cti  —  O  •  -  5  Av  |  —  4 .  i  - - 7* 


€ 

= 


m- 


ni  ^  ^ „  r- 

”  i»  c  ^ 


J  meujce 


p~=  ~y{Wxx  +  Wyy)  "  T°Vxx  ~ 

Pyy—  ~Y{Wxx  +  W™)  +  T(^XX  “  Wyv)  +  T?3^xy; 

pry=  pyx  =  -r?1VVry  +  -  Wyy); 


PXz  —  Pzx  =  —  TJoWxs  —  Tj^\Vyz\ 
Pyx  =  Pry  =  “^Wyr  +  V^xs', 
P  z  z  ■—  -VoWzz 

(here  the  z  axis  is  defined  parallel  to  B); 

SnfcT, 


ion  viscosity 


t7q  =  0.96n&TiTi;  = 


i  nfcT, 

73  _ 


2  Wei 


74  = 


nkTi 


IOuJtP 


V2  = 


6nkTi 


w. 


electron  viscosity 


..  nkTe  .  nkTe 

Vo  =  0.73nfcrere;  t;*=0.51  — 7a  =  2.0 


nfcT. 


73  =  “ 


.74  =  ~ 


nkTe 


“clT* 


7 

a;  J  r- 
ce  e 


2wce  wc_ 

For  both  species  the  rate-of-strain  tensor  is  defined  as 


dvj  dvk  2 

^=a^+a7-S^v'v- 


When  B  =  0  the  following  simplifications  occur: 


Ru  =  Tie}/ a  || ;  Rr  =  -0.71  nV(kTt);  q,  =  -*j|  V(A:T,); 

q'  =  0.71ni:reu:  q^-  =  -«.Jj  V{kTe );  Pjk  =  -T]0Wjk. 

For  wcere  ^  1  ^  wCI r, ,  the  electrons  obey  the  high-field  expressions  and  the 
ions  obey  the  zero-field  expressions. 

Collisionai  transport  theory  is  applicable  when  (1)  macroscopic  time  rates 
of  change  satisfy  djdt  <  1  /r,  where  r  is  the  longest  collisionai  time  scale,  and 
(in  the  absence  of  a  magnetic  field)  (2)  macroscopic  length  scales  L  satisfy  L 
l ,  where  l  =  vt  is  the  mean  free  path.  In  a  strong  .field.  wcer  ;>  1.  condition 
(2)  is  replaced  by  L y  l  and  £»x  VTr^  (Lx  >  r,  iu  a  uniform  field). 
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where  L,|  is  a  macroscopic  scale  parallel  to  the  field  B  and  L j_  is  the  smaller 
of  B/\V±B\  and  the  transverse  plasma  dimension.  In  addition,  the  standard 
transport  coefficients  are  valid  only  when  (3)  the  Coulomb  logarithm  satisfies 
A  >  1:  (4)  the  electron  gyroradius  satisfies  re  >  \D ,  or  3rrnemec 2  >  B2;  (5) 
relative  drifts  u  =  va  -  vj  between  two  species  are  small  compared  with  the 
thermal  velocities,  i.e.,  u  kTa  / ma ,  kTj/vnj\  and  (6)  anomalous  transport 
processes  owing  to  microinstabilities  are  negligible. 

Weakly  Ionized  P  las  mas 

Collision  frequency  for  scattering  of  charged  particles  of  species  c*  by 
neutrals  is 

*<*  =  n0a°/0(kTa/ma)1/2, 

where  n0  is  the  neutral  density  and  is  the  cross  section,  typically  ^ 

5  x  10  15  cm"  and  weakly  dependent  on  temperature. 

When  the  system  is  small  compared  with  a  Debye  length,  L  <  \D,  the 
charged  particle  diffusion  coefficients  are 


Da  =  kTa  /mQ  i/a  , 

In  the  opposite  limit,  both  species  diffuse  at  the  ambipolar  rate 

D  =  fiiDe  ~  =  (Ti  +  T')D'D e 

M»  Me  T,  De  +  TeDi 

where  =  ea /m0  i/a  is  the  mobility.  The  conductivity  crQ  satisfies  <rQ  = 

^  ry  €  C*  f-L  o  • 

In  the  presence  of  a  magnetic  field  B  the  scalars  /j,  and  <r  become  tensors, 
J*  =  a*  ■  E  =  E.J  +  j"Ex  +  <r®  E  x  b. 

where  b  =  B/B  and 

<7^  =  Tie,  Co2 /ma  v<y  - 
ct  a  2  , ,  2  2  . 

17  ±-  _  <7;|  Vr*  /(^or  +WCJ; 

ay\  —  ^||  *JCa  /(  Vf,  f  U/CQ  )  . 

Here  a  ±  and  erA  are  the  Pedersen  and  Hall  conductivities,  respectively. 


39 


APPROXIMATE  MAGNITUDES 
IN'  SOME  TYPICAL  PLASMAS 


Plasma  Type 

_  3 

a  cm 

T  eV 

jpe  sec  1 

A  &  c:n 

n  Xd2 

vei  sec”  k 

Interstellar  gas 

1 

1 

6  x  104 

7  x  102 

4  x  10s 

7  x  10-5 

Gaseous  nebula 

103 

1 

2  x  106 

20 

107 

6  x  10-2 

Solar  Corona 

106 

102 

6  x  10~ 

7 

4  x  103 

6  x  10-2 

Diffuse  hot  plasma 

1012 

102 

6  x  1010 

7  x  10-3 

4  x  105 

40 

Solar  atmosphere. 

1014 

1 

6  x  10u 

7  x  10-s 

40 

2  x  109 

gas  discharge 

Warm  plasma 

1014 

10 

6  x  10X1 

2  x  10-4 

103 

107 

Hot  plasma 

10 14 

102 

6  x  1011 

7  x  10~4 

4  x  104 

4  x  106 

Thermonuclear 

1015 

104 

2  x  1012 

2  x  10"3 

10T 

5  x  104 

plasma 

Theta  pinch 

1016 

102 

6  x  1012 

7  x  10_s 

4  x  103 

3  x  10s 

Dense  hot  plasma 

1013 

102 

6  x  1013 

7  x  10~c 

4  x  102 

2  x  1010 

Laser  Plasma 

h-* 

o 

IJ 

o 

102 

6  x  1014 

7  x  10_r 

40 

2  x  1012 

The  diagram  (facing)  gives  comparable  information  in  graphical  form.22 


40 


-2-10  I  2  3  4  5 

Ioqiq  T(eV) 


41 


IONOSPHERIC  PARAMETERS22 


The  following  tables  give  average  nighttime  values.  Where  two  num¬ 
bers  are  entered,  the  first  refers  to  the  lower  and  the  second  to  the  upper 
portion  of  the  layer. 


Quantity 

E  Region 

F  Region 

Altitude  (km) 

90-160 

160-500 

Number  density  (m-3) 

1.5  x  1010-3.0  x  1010 

5  x  1010-2  x  1011 

Height-integrated  number 
density  (m-2) 

9  x  1014 

4.5  x  1015 

Ion-neutral  collision 
frequency  (sec-1) 

2  x  103-102 

0.5-0.05 

Ion  gyro-/colIision 
frequency  ratio  ac  , 

0.09-2.0 

4.6  x  102-5.0  x  103 

Ion  Pederson  factor 

Ki/(1  +  AC  i2) 

0.09-0.5 

2.2  x  10-3-2  x  10-4 

Ion  Hall  factor 

8  x  10-4-0.8 

1.0 

Electron-neutral  collision 
frequency 

1.5  x  104-9.0  x  102 

80-10 

Electron  gyro-/collision 
frequency  ratio  ac* 

4.1  X  102-6.9  x  103 

7.8  x  104-6.2  x  10s 

Electron  Pedersen  factor 

ACe/(l  +  At*2) 

2.7  x  10-3-r.5  x  10-4 

10-s-1.5  x  10-6 

Electron  Hall  factor 

ac*2/(1  +  AC*2) 

1.0 

1.0 

Mean  molecular  weight 

28-26 

22-16 

Ion  gyrofrequency  (sec-1) 

ISO-190 

230-300 

Neutral  diffusion 

coefficient  (m"  sec-1) 

30-5  x  103 

10s 

The  terrestrial  magnetic  field  in  the  lower  ionosphere  at  equatorial 
lattitudes  is  approximately  Bo  =  0.35  x  10  1  tesla.  The  earth’s  radius  is 

Re  =  6371  km. 
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SOLAR  PHYSICS  PARAMETERS24 


Chromosphere  and  Corona25 


Parameter  (Units) 

Quiet 

Sun 

Coronal 

Hole 

Active 

Region 

Chromospheric  radiation  losses 
(erg cm-2  s~ 1 ) 

Low  chromosphere 

2  x  106 

2  x  106 

£  io7 

Middle  chromosphere 

2  x  106 

2  x  106 

107 

Upper  chromosphere 

3  x  105 

3  x  105 

2  x  10G 

Total 

4  x  106 

4  x  106 

>  2  x  107 

Transition  layer  pressure  (dyne  cm-2) 

0.2 

0.07 

2 

Coronal  temperature  (K)  at  1.1  R,t> 

1.1-1. 6  x  10s 

10G 

2.5  x  10G 

Coronal  energy  losses  (erg  cm-2  s-1) 

Conduction 

2  x  10s 

G  x  104 

105-107 

Radiation 

10s 

104 

5  x  10G 

Solar  Wind 

<  5  x  104 

7  x  10s 

<  10s 

Total 

3  x  105 

8  x  10s 

107 

Solar  wind  mass  loss  (g  cm-2  s-1 ) 

<  2  x  10“11 

2  x  10-l° 

<  4  x  10“11 
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THERMONUCLEAR  FUSION20 

Natural  abundance  of  isotopes: 

hydrogen  n D/nH  =  1.5  x  10~4 
helium  nHe3  /nHe4  =  1-3  x  10“° 

lithium  nLi6/nLi7  =  0-0S 

Mass  ratios:  meJm.D  —  2.72  x  10~4  =  1/3G70 

(me/mD)1/2  =  1.65  X  ICC2  =  1/60.6 
m(/mr  =  1.82  X  10~4  =  1/5496 
(me/mT)1/2=  1.35  X  10~2  =  1/74.1 

Absorbed  radiation  dose  is  measured  in  rads:  1  rad  =  102ergg_l.  The 
curie  (abbreviated  Ci)  is  a  measure  of  radioactivity:  1  curie  =  3.7  x 
1010  counts  sec-1. 

Fusion  reactions  (branching  ratios  are  correct  for  energies  near  the  cross 
section  peaks;  a  negative  yield  means  the  reaction  is  endothermic):2' 


(la) 

D  +  D  50% 

►T(l.OlMeV)  +  p(3.02  MeV) 

(lb) 

50% 

>He3(0.82  MeV)  +  n(2.45MeV) 

(2) 

D  +  T  - 

►He4(3.5MeV)  +  u(14.1MeV) 

(3) 

D  +  He3 - 

►He4(3.6  MeV)  +  p(14.7MeV) 

(4) 

T  +  T  - 

►He4  +  2u  +  11.3  MeV 

(5a) 

HeI  +  T  51% 

►He4  +  p  +  n  +  12.1  MeV 

(5b) 

43% 

►He4(4.8  MeV)  +  D(9.5  MeV) 

(5c) 

6% 

►He5(2.4  MeV)  +  p(11.9MeV) 

(6) 

p  +  Li6  - 

►He4(1.7  MeV)  +  He3(2.3MeV) 

(7a) 

P  +  Ll  20% 

♦2  He4  -f  17.3  MeV 

(7b) 

80% 

♦Be'  +  u  —  1.6  MeV 

(8) 

D  +  Li6  - 

*2He4  4-  22.4  MeV 

O) 

p  +  B11  - 

■*3  He4  +  8.7  MeV 

(10) 

n  +  Li6 

*He4(2.1  MeV)  +  T(2.7MeV) 

The  total  cross  section  in  barns  as  a  function  of  E,  the  energy  in  kcV 
of  the  incident  particle  [the  first  ion  on  the  left  side  of  Eqs.  (l)-(5)], 
assuming  the  target  ion  at  rest,  can  be  fitted  by23 

As,  +  f(A4  -  AzE)2  +  l]  A2 

*t(E)  =  - V - ^ - 

E  exp(AlE-V 2)  -  1 
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where  the  Duane  coefficients  Aj  for  the  principle  fusion  reactions  are  as  follows: 


1 

D-D 

(la) 

D-D 

(lb) 

D-T 

(2) 

D-He3 

(3) 

T-He3 

(5a-c) 

Ax 

46.097 

47.88 

45.95 

89.27 

38.39 

123.1 

372 

482 

50200 

25900 

44S 

11250 

^3 

4.36  x  10-4 

3.08  x  10-4 

1.368  x  10-2 

3.98  x  10-3 

1.02  x  10-3 

0 

1.220 

1.177 

1.076 

1.297 

2.09 

0 

As 

0 

0 

409 

647 

0 

0 

Reaction  rates  av  (in  cm3  sec-1),  averaged  over  Maxwellian  distributions: 


Temperature 

(keV) 

D-D 

(la  +  lb) 

D-T 

(2) 

D-He3 

(3) 

wm 

T-He3 

(5a-c) 

1.0 

1.5  x  10-22 

5.5  x  10-21 

10  — 28 

3.3  x  10-22 

IQ-28 

2.0 

5.4  x  10-21 

2.6  x  10-19 

1.4  x  10-23 

7.1  x  10-21 

10-25 

5.0 

1.8  x  10-19 

1.3  x  10“17 

6.7  x  10-21 

1.4  x  10-19 

2.1  x  10-22 

10.0 

1.2  x  10-18 

1.1  x  10-16 

2.3  x  10-19 

7.2  x  10-19 

1.2  x  10-2° 

20.0 

5.2  x  10~18 

4.2  x  10-16 

3.8  x  10-18 

2.5  x  10-18 

2.6  x  10-19 

50.0 

2.1  x  10-17 

8.7  x  10-16 

5.4  x  10-17 

8.7  x  10-18 

5.3  x  10-18 

100.0 

4.5  x  10-17 

8.5  x  10-16 

1.6  x  10-16 

1.9  x  10-17 

2.7  x  10-17 

200.0 

8.8  x  10-17 

6.3  x  10-18 

2.4  x  10-16 

4.2  x  10-17 

9.2  x  10-17 

500.0 

1.8  x  10-16 

3.7  x  10-16 

2.3  x  10-16 

8.4  x  10-17 

2.9  x  10-16 

1000.0 

2.2  x  10-16 

2.7  x  10-16 

1.8  x  10-16 

8.0  x  10-17 

5.2  x  10-16 

For  low  energies  (T  ^  25  keV)  the  data  may  be  represented  by 

Rdo  =  2.33  x  10”l4r_2/3  exp(-18.76T-1/3)  cm3  sec-1; 


( <rv)DT  =  3.68  x  10-12T-2/3  exp(-19.94T-1/3)  cm3  sec-1, 
where  T  is  measured  in  keV. 

The  power  density  released  in  the  form  of  charged  particles  is 

Pdc  —  3.3  x  10-l3nD2(  w)dd  watt  cm  3  (including  the  subsequent 

D-T  reaction); 

Pdt  =  5.6  x  10-13 nonT^Pv) dt  watt  cm-3; 

PDH*Z  —  2-9  *  10"12jlDnHe3(^)DHe3  wattcm-3. 
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RELATIVISTIC  ELECTRON  BEAMS 


Here  7  =  (1  —  /?2)-1//2  is  the  relativistic  scaling  factor;  quantities 
in  analytic  formulas  are  expressed  in  SI  or  cgs  units,  as  indicated;  in 
numerical  formulas.  I  is  in  amperes  (A),  S  is  in  gauss  (G),  electron 
linear  density  .V  is  in  cm-1,  and  temperature,  voltage  and  energy  are  in 
MeV;  0Z  =  vx/c\  k  is  Boltzmann’s  constant. 

Relativistic  electron  gyroradius: 


mc~ 

~7b 


(7"  —  l)1^'  (cgs)  =  1.70  x  10J(7*  —  l)l,~  B~l  cm 


3/2  _  i\l/2  n-1 


Relativistic  electron  energy: 


W  =  mc2-j  =  0.5117  MeV. 


Bennett  pinch  condition; 


I2  =  2Nk(Te  +  T,)c 2  (cgs)  =  3.20  x  10"4.V(Te  +  T, )  A2. 
Alfven-Lawson  limit: 

IA  =  (mc3/e)/?r 7  (cgs)  =  (4irmc/ noe)0x-r  (SI)  =  1.70  x  104/3,7A. 
The  ratio  of  net  currer. .  to  IA  is 


J_  _  u_ 

I A  7  ’ 

Here  1/  =  iVre  is  the  Budker  number,  where  re  =  e2/mc2  =  2.82  x 
10" 13  cm  is  the  classical  electron  radius.  Beam  electron  number  density 
is 

nb  =  2.08  x  1087^?-1  cm-3, 

where  J  is  the  current  density  in  A  cm-2.  For  a  uniform  beam  of  radius 
a  (in  cm), 

rib  =6.63  x  10'  Ia~2 0~l  cm-3 , 

and 

2re  1/ 

a  7  ‘ 
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Child’s  law:  (non*relativistic)  space-charge-limited  current  density  be¬ 
tween  parallel  plates  with  voltage  drop  V  and  separation  d  (in  cm)  is 

J  -  2.34  x  10 3V3/2d~2  A  cm-2. 

The  saturated  parapotential  current  (magnetically  self-limited  flow  along 
equipotentials  in  pinched  diodes  and  transmission  lines)  is20 

Jp  =  8.5  x  103C?7  In  [7  +  (72  -  1)1/2]  A. 

where  G  is  a  geometrical  factor  depending  on  the  diode  structure: 

for  parallel  plane  cathode  and  anode 
of  width  w ,  separation  d; 

for  cylinders  of  radii  R\  (inner)  and  H2  (outer); 

for  conical  cathode  of  radius  Rc ,  maximum 
separation  do  (at  r  =  Rc)  from  plane  anode. 

For  (3  — *  0  (7  — ►  1),  both  1a  and  lp  vanish. 

The  condition  for  a  longitudinal  magnetic  field  B-  to  suppress  filamen- 
tation  in  a  beam  of  current  density  J  (in  A  cm-2)  is 

B,  >  47/3r(7J)1/2G. 

Voltage  registered  by  Rogowski  coil  of  minor  cross-sectional  area  A ,  n 
turns,  major  radius  a ,  inductance  L ,  external  resistance  R  and  capaci¬ 
tance  C  (all  in  SI): 

externally  integrated  V  =  ( 1/ RC)(nA/io  //2xa); 

self-integrating  V  =  (R/ L){nApo  I /2na)  =  Rijn. 

X-ray  production,  target  with  average  atomic  number  Z  ( V  ^  5  MeV): 
tj  =  x-ray  power/beam  power  =  7  x  10 -4ZV\ 

X-ray  dose  at  1  meter  generated  by  an  e-beam  depositing  total  charge  Q 
coulombs  while  V  >  0.84V'm*x  in  material  with  charge  state  Z : 

D  =  150  Vm2x8  QZ1/2  rads. 
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BEAM  INSTABILITIES30 


Name 

Conditions 

Saturation  Mechanism 

Electron- 

electron 

Vi  >  Vej,  j  =  1.2 

Electron  trapping  until 

vej  ~  Vd 

Buneman 

Vi  >  (A//m)1/3Vi, 
V*  >  Ve 

Electron  trapping  until 

Ve  -  Vd 

Beam-plasma 

vb  >  (np/nb)l/3Vb 

Trapping  of  beam  electrons 

Weak  beam- 
plasma 

Vb  <  (np/nb)1/3Vb 

Quasiliuear  or  nonlinear 
(mode  coupling) 

Beam-plasma 

(hot-electron) 

ve>vb>  Vb 

Quasilinear  or  nonlinear 

Ion  acoustic 

Te  >  Ti,  Vd  >  C, 

Quasiliuear,  ion  tail  form¬ 
ation,  nonlinear  scattering, 
or  resonance  broadening. 

Anisotropic 

temperature 

(hydro) 

TeX  >  2T.ii 

Isotropization 

Ion  cyclotron 

Vd  >  20 Vi  (for 

T'VTi) 

Ion  heating 

Beam-cyclotron 

(hydro) 

vd  >  C. 

Resonance  broadening 

Modified  two- 
stream  (hydro) 

Vd  <(1  -f /?)l/2VA, 
Vi  >  C. 

Trapping 

Ion-ion  (equal 
beams) 

U  <  2(1  +  f3)l/2VA 

Ion  trapping 

Ion-ion  (equal 
beams) 

U  <  2C, 

Ion  trapping 

For  nomenclature,  see  p.  50. 
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Parameters  of  Most  Unstable  Mode 


Name 

Wave 

Group 

Growth  Rate 

Frequency 

Number 

Velocity 

1 

we 

Electron- 

0 

0.9-^ 

0 

electron 

2 

Vd 

(m  Y'3 

f  m  \  1/3 

we 

2 

Buneman 

0.7  —  ue 

0.4  [  —  )  we 

-Vd 

\Mj 

Vd 

3 

Beam-plasma 

0.7  f  — )  '  we 

we  - 

We 

|vb 

\np) 

Vb 

3 

Weak  beam- 

nb  (Vb\2 

2 np  VvJ  We 

<jJe 

we 

3  Ve2 

plasma 

Vb 

n 

Beam-plasma 

(liot-electron) 

(  \  V, 

UJ 

Vb 

—  ue 

Ve 

AoJ 

Vb 

Ion  acoustic 

(5)"- 

Wi 

c. 

Anisotropic 

we  cos  0  ~  f2e 

-r1 

Vex 

temperature 

(hydro) 

Ion  cyclotron 

0.1  Qi 

1.2ft,- 

r~l 

1 

5* 

Beam-cyclotron 

o.7ne 

nfte 

0.7A-1 

£  Vd; 

(hydro) 

nH 

1.7  — 

<  c. 

Modified  two- 

“ &H 

0.9  ftH 

stream 

2 

Vd 

2 

(hydro) 

n„ 

1.2  — 

Ion-ion  (equal 

0.4Qh 

0 

0 

beams) 

£/ 

Ion-ion  (equal 

0.4w; 

0 

U>i 

l-217 

0 

beams) 

u 

For  nomenclature,  see  p.  50. 
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In  the  preceding  tables,  subscripts  e.  t.  d.  b.  p  stand  for  "electron.” 
“ion.”  "drift.”  "beam. ”  and  "plasma.”  respectively.  Thermal  velocities 
are  denoted  by  a  bar.  In  addition,  the  following  are  used: 


m  electron  mass  re . 

A f  ion  mass  J 

V  velocity 

T  temperature  V± 

ne,n,  number  density  f2e.Q, 

n  harmonic  number  Qh 

C,  —  (Te  / ion  sound  speed 
we,  w,  plasma  frequency  U 

A d  Debye  length 


gvroradius 

plasma/ magnetic  energy 
density  ratio 
Alfveu  speed 
gyrofrequeucy 
hybrid  gyrofrequeucy. 

nH-  =  QeQ, 
relative  drift  velocity  of 
two  ion  species 


LASERS 


System  Parameters 


Efficiencies  and  power  levels  are  approximately  state-of-the-art  (1987). 31 


Type 

Wavelength 

Efficiency 

Power  levels  available  (W) 

(pm) 

Pulsed 

cw 

O 

O 

10.6 

0.01-0.02 

(pulsed) 

>  2  x  1013 

>  105 

CO 

5 

0.4 

>  109 

>  100 

Holmium 

2.06 

0.03 

>  107 

30 

Iodine 

1.315 

0.003 

>  1012 

- 

Nd-glass, 

YAG 

1.06 

0.001-0.06 

~  1014 

(10-beam  system) 

1-103 

*Color 

center 

1-4 

10-3 

>  106 

1 

*OPO 

0. 7-0.9 

10“3 

106 

1 

Ruby 

0.6943 

<  10~3 

O 

o 

H 

1 

He-Ne 

0.6328 

10-4 

- 

1-50  x 10-3 

*Argon  ion 

0.45-0.60 

10~3 

5  x  104 

1-10 

n2 

0.3371 

0.001-0.05 

105-106 

- 

*Dye 

0.3-1. 1 

10-3 

>  106 

140 

Kr-F 

0.26 

0.08 

>  109 

- 

Xenon 

0.175 

0.02 

>  10s 

- 

^Tunable  sources 

YAG  stands  for  Yttrium-Aluminum  Garnet  and  OPO  for  Optical  Para¬ 
metric  Oscillator. 
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Formulas 

An  e-m  wave  with  k  |]  B  has  an  index  of  refraction  given  by 
n±  =  [1  -  upl/u(u  ^  wce)]1/2, 

where  ±  refers  to  the  helicity.  The  rate  of  change  of  polarization  angle 
9  as  a  function  of  displacement  s  (Faraday  rotation)  is  given  by 

dd/ds  =  (fc/2)(n_  -  n+)  =  2.36  x  10 4NBf~2  cm"1, 

where  N  is  the  electron  number  density,  B  is  the  field  strength,  and  /  is 
the  wave  frequency,  all  in  cgs. 

The  quiver  velocity  of  an  electron  in  an  e-m  field  of  angular  fre¬ 
quency  u  is 


vo  =  eEm±x /rnu>  —  25.6I1/,2Ao  cm  sec"1 

in  terms  of  the  laser  flux  I  =  cFm2x/87r,  with  I  in  watt/cm2,  laser 
wavelength  A0  in  y.m.  The  ratio  of  quiver  energy  to  thermal  energy  is 

W^/Wth  =  rnev07/2kT  -  1.81  x  10"13Ao2//T, 

where  T  is  given  in  eV.  For  example,  if  /  =  1015  Wcm'2,  A0  =  1  /im, 
T  =  2keV,  then  WqU/Wth  »  0.1. 

Pondermotive  force: 

T  -  NV{E2)/%*Nc, 

where 

Nc  -  1.1  x  1021  Ao  "2cm”3. 

For  uniform  illumination  of  a  lens  with  /-number  F,  the  diameter 
d  at  focus  (85%  of  the  energy)  and  the  depth  of  focus  l  (distance  to  first 
zero  in  intensity)  are  given  by 

d  «  2A4F\9/6dl  and  l  a  ±2F2 \9/9DL . 

Here  9  is  the  beam  divergence  containing  85%  of  energy  and  9ql  is  the 
diffraction-limited  divergence: 

9dl  =  2.44A/6, 

where  b  is  the  aperture.  These  formulas  are  modified  for  nonuniform 
(such  as  Gaussian)  illumination  of  the  lens  or  for  pathological  laser  pro¬ 
files. 
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ATOMIC  PHVSICS  AND  RADIATION 


Energies  and  temperatures  are  in  eV:  all  other  units  are  cgs  except 
where  noted.  Z  is  the  charge  state  (Z  =  0  refers  to  a  neutral  atom); 
the  subscript  e  labels  electrons.  N  refers  to  number  density,  n  to  princi¬ 
pal  quantum  number.  Asterisk  superscripts  on  level  population  densities 
denote  local  thermodynamic  equilibrium  (LTE)  values.  Thus  Nnm  is  the 
LTE  number  density  of  atoms  (or  ions)  in  level  n. 

Characteristic  atomic  collision  cross  section: 

(1)  ira02  =  8.80  x  10-17  cm2. 


Binding  energy  of  outer  electron  in  level  labelled  by  quantum  numbers 
n,  l : 


(2) 


£»(n,/)  = 


Z2Et 


(n  —  Aj)2  ’ 


where  E ^  =  13.6 eV  is  the  hydrogen  ionization  energy  and  A j  =  0.751-5 
1  >  5,  is  the  quantum  defect. 


Excitation  and  Decay 

Cross  section  (Bethe  approximation)  for  electron  excitation  by  dipole 
allowed  transition  m  — *  n  (Refs.  32,  33): 


(3) 


^mn 


=  2.36  x  10~13 


fnmQ^Tli  m) 
«A  Erxm 


cm", 


where  /nm  is  the  oscillator  strength,  g(n ,  m)  is  the  Gaunt  factor,  e  is  the 
incident  electron  energy,  and  A Enm  ~  En  —  Em- 

Electron  excitation  rate  averaged  over  Maxwellian  velocity  distribution, 
Xmn  =  Nt(<7mnv)  (Refs.  34,  35): 


(4) 


Xmn  =  1-6  X  10 


—  5  f  nm  ( g(n ,  m))Ne 


A£nm77 


1/2 


exp 


A£nm  \ 

J 


sec  , 


where  (g(n,  m))  denotes  the  thermal  averaged  Gaunt  factor  (generally  ~  1 
for  atoms,  ~  0.2  for  ions). 
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Rate  for  electron  collisional  deexcitation: 


(5)  ynm  =  (Nm  m/Nn*)XTnn. 

Here  Nm*/Nn*  =  {Om  /  On)  expfAEnm  /Te)  is  the  Boltzmann  relation  for 
level  population  densities,  where  gn  is  the  statistical  weight  of  level  a. 

Rate  for  spontaneous  decay  n  —  m  (Einstein  A  coefficient)34 

(6)  inm  =  4.3  X  10‘  (gn  /gm)fnm(&£nm)2  Sec"1. 


Intensity  emitted  per  unit  volume  from  the  transition  n  — •  m  in  an  optic¬ 
ally  thin  plasma: 

(7)  Inm  —  1.6  x  10 AnmNn^Erxm  watt/cm3. 


Condition  for  steady  state  in  a  corona  model: 

(8)  NoNe(<ronv)  =  NnAn0, 


where  the  ground  state  is  labelled  by  a  zero  subscript. 
Hence  for  a  transition  n  — »  m  in  ious,  where  (g(n,  0))  s; 


(9) 


=  5.1  x  10 


—  25  fnmgoNeNo 


9  m  Tg 


1/2 


0.2. 

A£n  o\ 

T*  ) 


watt 

cm3 


Ionization  and  Recombination 

In  a  general  time-dependent  situation  the  number  density  of  the 
charge  state  Z  satisfies 


(10) 


dN(Z) 

dt 


Ne 


-  S(Z)N(Z)  -  a(Z)N(Z) 


+  S{Z  -  1  )N{Z  -  1)  +  oc(Z  4-  1  )N(Z  4-  1) 


Here  S(Z)  is  the  ionization  rate.  The  recombination  rate  a(Z)  has  the 
form  oc(Z)  —  ar(Z)  4-  Neaj(Z),  where  ar  and  az  are  the  radiative  and 
three-body  recombination  rates,  respectively.* 
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Classical  ionization  cross-section36  for  any  atomic  shell  j 


'll' 


<r,  =  6  x  10  14  bj  g:  (z)/lT,  ~  cm' . 


Here  b:  is  the  number  of  shell  electrons:  U :  is  the  binding  energy  of  the 
ejected  electron:  x  =  e/U where  e  is  the  incident  electron  energy:  and  g 
is  a  universal  function  with  a  minimum  value  gmia  S:  0.2  at  z  %  4. 

Ionization  from  ion  ground  state,  averaged  over  Maxwellian  electron  dis¬ 
tribution.  for  0.02  <  Te,‘  <  100  (Ref.  35): 


(12) 


S(Z)  =  10‘ 


(Te/E*)1/- 


(£^)3/-(o.oTre/f£ 


exp  - 


Te 


3  , 

cm  / sec. 


where  Et^  is  the  ionization  energy. 

Electron-ion  radiative  recombination  rate  (e  -f-  .V(Z)  —  N(Z  —  1)  —  h.v) 
for  Te/Z2  <  400  eV  (Ref.  37): 


(13)  ar(Z)  =  5.2  x  10 


Ez  N  1/^2 

-14Z  f 

Te 


0.43-  -In(E^/Te) 


+0.469(  E^/Te )  ~ 1/3 


3  / 

cm  /sec. 


For  1  eV  <  Te/Z2  <  15  eV,  this  becomes  approximately35 

(14)  ar{Z)  =  2.7  x  10“13Z2Te"l/2  cm3/sec. 


Collisional  (three-body)  recombination  rate  for  singly  ionized  plasma:38 

(15)  a3=S.75xl0-*'Te~4'5cm6/sec. 


Photoionization  cross  section  for  ions  in  level  n.  I  (short-wavelength  limit): 

(16)  <Tph( n.  1)  =  1.64  x  10-16Z°/n3K"'!"‘;cm‘, 

where  K  is  the  wavenumber  in  Rydbergs  (1  Rydberg  =  1.0974  x  10°  cm-1 ). 
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Ionization  Equilibrium  Models 
Saha  equilibrium:39 


(17) 


Ne  N\*(Z) 
Nn*(Z  -  1) 


=  6.0  x  10' 


:i  Q\  Te3^2 


.2  —  1 
9  n 


exp  - 


Ej(n.l) 

Te 


cm 


-3 


where  is  the  statistical  weight  for  level  n  of  charge  state  Z  and  E^^n.  1) 
is  the  ionization  energy  of  the  neutral  atom  initially  in  level  (n,  i),  given 
by  Eq.  (2). 

In  a  steady  state  at  high  electron  density, 


(18)  NeN*(Z)  S(Z  —  1) 

1  ’  N*(Z  -  1)  a3 

a  function  only  of  T. 

Conditions  for  LTE:39 

(a)  Collisioaa.1  and  radiative  excitation  rates  for  a  level  n  must  satisfy 

(19)  Vnm  £  lOAnm* 


(b)  Electron  density  must  satisfy 

(20)  N€  >  7  x  10l8Z7n”l7/2(T/£*  )1/2cm-3. 


Steady  state  condition  in  corona  model: 


(21) 


N{Z  -  1)  Qr 

N{Z)  "  S(Z  -  1)' 


Corona  model  is  applicable  if40 

(22)  1012t/_1  <  Ne  <  1016T«7/2  cm-3, 

where  tj  is  the  ionizatiou  time. 
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Radiation 

N.  B.  Energies  and  temperatures  are  in  eV:  all  other  quantities  are 
in  cgs  units  except  where  noted.  Z  is  the  charge  state  i  Z  =  0  refers  to  a 
neutral  atom);  the  subscript  e  labels  electrons.  .V  is  number  density. 

Average  radiative  decay  rate  of  a  state  with  principal  quantum  number  u 
is 


(23) 


=  l.b  x  10  Z  n 


m  <  n 


sec. 


Natural  linewidth  (A E  in  eV): 

(24)  A£  At  =  h  =  4.14  x  10” 15  eV  sec. 

where  A f  is  the  lifetime  of  the  line. 

Doppler  width: 

(25)  AA/A  =  7.7xlO"5(r//r)l/:, 

where  /i  is  the  mass  of  the  emitting  atom  or  ion  scaled  by  the  proton  mass. 
Optical  depth  for  a  Doppler-broadened  line:39 

(26)  r  =  1.76  x  10" 13  c2 / kT)1/2  NL  -  5.4  x  10~9  \((i/T)l/2  N  L. 

where  A  is  the  wavelength  and  L  is  the  physical  depth  of  the  plasma;  M, 
JV,  and  T  are  the  mass,  number  density,  and  temperature  of  the  absorber; 
/x  is  M  divided  by  the  proton  mass.  Optically  thin  means  r  <  1. 

Resonance  absorption  cross  section  at  center  of  line: 

(27)  =  5.6  x  10-13 A2/AA  cm2. 

Wien  displacement  law  (wavelength  of  maximum  black-body  emission): 
(2S)  Amax  =  2.50  x  10-5T-1  cm. 

Radiation  from  the  surface  of  a  black  body  at  temperature  T: 

(29)  W  =  1.03  x  105T4  watt/cm2. 
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Bremsstrahlung  from  hydrogen-like  plasma:26 

(30)  Pbt  =  1-69  x  10"32  iVeTe1/2  ^  [z2iV(Z)]  watt/ 


cm 


where  the  sum  is  over  all  ionization  states  Z. 

Bremsstrahlung  optical  depth:41 

(31)  r  =  5.0  x  lQ-38NcNiZ2gLT~7/2, 

where  ~g  ss  1.2  is  an  average  Gaunt  factor  and  L  is  the  physical  path  length. 

Inverse  bremsstrahlung  absorption  coefficient42  for  radiation  of  angular 
frequency  u: 

(32)  «  =  3.1X  10“7Zne2  lnAT"3/V2(l  -  wj/u,2)1/2  cm-1; 

here  A  is  the  electron  thermal  velocity  divided  by  V,  where  V  is  the  larger 
of  u)  and  up  multiplied  by  the  larger  of  Ze?  /kT  and  K/{mkT)xl2 . 

Recombination  (free-bound)  radiation: 


(33)  Pr  =  1.69  x  10'32^ere1/2  ^  Z2N{Z)  watt/cm3 


Cyclotron  radiation28  in  magnetic  field  B: 

(34)  Pe  -  6.21  X  10~28 J9 2 IVeTe  watt/cm3. 


For  NekTe  =  NikTi  =  B2 /16*  (/?  =  1,  isothermal  plasma),26 

(35)  Pc  =  5.00  x  10 _38N2T2  watt/cm3. 


Cyclotron  radiation  energy  loss  e-folding  time  for  a  single  electron:41 
(36) 


9.0  x  10aB"2 
tc  «  - — - — ; - sec. 


2.5  +  7 


where  7  is  the  kinetic  plus  rest  energy  divided  by  the  rest  energy  me2. 
Number  of  cyclotron  harmonics41  trapped  in  a  medium  of  finite  depth  L : 

(37)  mtr  =  (570BL)1/6, 

where  0  =  SirNkT/ B2 . 

Line  radiation  is  given  by  summing  Eq.  (9)  over  all  species  in  the  plasma. 
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